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Numerical  results  are  present eJ  for  the  approximate 
ime-averaged  power  density  dissipated  on  the  surface  of 
highly  conductive  infinite  wedge  (of  angle  90°)  near 
he  vicinity  of  the  sharp  edge.  I  lie  numerical  results 
re  normalized  so  one  may  readily  obtain  results  for 
trying  incident  wave  frequencies.  The  angle  of  incidence 
:  the  incident  radiation,  0(  (measured  with  respect  to 
,e  surface  normal  of  the  wedge),  is  varied  from  Q*  -  0° 
normal  incidence)  to  =  90°  (grazing  incidence).  In 
irticular,  numerical  results  for  IT  and  TM  polarization 
re  presented  for  incident  wavelengths  of  3 . 8  uni  and 
r.b^un  .  It  is  assumed  that  the  power  density  absorbed 

y 

at!  /in")  is  all  converted  to  heat  density. 

iv  accept  the  validity  of  I  lie  results,  an  extensive 
ec is  made  on  ,ui  arbitrary  (real)  order  Bessel  function 
■threat  ine  used  in  the  equation  necessa  r for  the  numerical 
.i  1  c>.  1  at  i  me  .  tieneral  e  lect  ro»<ngnet  i  c  'iieory  is  also 
iscu.-sed  and  extended  to  the  use  in  specific  geometries, 
fe  geometries  considered  besides  the  wedge  (though  with 
e  numerical  results)  are  the  parabolic  cylinder  and 
1 1  f in : t e  cone. 

In  conclusion,  the  numerical  results  seem  very  promising 
,.e  to  the  success  of  verifying  the  Bessel  function  sub- 
out  ine  used  in  the  calculations. 
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INVESTIGATION  OF  LASER  HEATING 
IN  COMPONENTS  WITH  SHARP  EDGES 

I .  Introduct ion 

The  number  of  lasers  being  used  has  rapidly  increased 
in  the  past  few  years  due  in  a  large  part  to  the  advent 
ol"  more  sophisticated  laser  systems  and  the  discovery  of 
new  laser  applications.  The  improvement  in  quality  of 
optical  components  in  lasers  has  also  increased.  As  a 
esult,  the  cost  of  many  of  the  high  quality  optical 
components  has  risen. 

Of  particular  interest  to  the  Air  Force  Weapons 
Laboratory  (AFWL)  at  Kirtland  Air  Force  Base,  New  Mexico, 
is  the  investigation  of  the  power  dissipated  per  unit 
area  on  the  surfaces  of  highly  conductive  dielectrics  in 
the  vicinity  of  sharp  edges  when  illuminated  by  high 
intensity  radiation.  The  sharp  edges  arise  from  the  use 
of  rectangular  mirrors  (having  sharp  edges)  and  of  cone- 
shaped  mirrors  used  as  waxicons  and  reflaxicons  (see 
Figures  3  and  4  in  Chapter  II).  Of  special  interest  to 
ARVL  is  optical  components  having  a  conductivity  near 
that  of  copper.  The  incident  radiation  of  interest  is 
in  the  micron  wavelength  range  (3.8yum  and  10.6^(m,  to 
he  exact).  It  is  feared  that  the  damage  threshold  of  these 
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opt ical  components  may  be  surpassed  near  the  vicinity  of 
the  sharp  edges  due  to  induced  surface  currents,  thus 
having  the  potential  of  damaging  expensive  optical  components. 

The  main  objective  of  this  thesis  is  to  calculate  the 
time-averaged  power  density  dissipated  on  the  surface  of 
a  highly  conductive  90°  wedge  imbedded  in  free  space.  It 
is  assumed  that  the  power  density  dissipated  represents 
the  heat  density  (in  watts  per  square  meter)  produced  in 
the  wedge  from  the  incident  field.  The  power  density  is 
found  by  first  assuming  the  wedge  to  be  a  perfect  conductor 
uno.  finding  the  surface  current  density  induced  on  the 
.av.,e  surface.  An  approximation  is  then  made  to  find  the 
power  density  using  the  skin  deptii  for  a  highly  conductive 
linear  medium  and  the  surface  current  density  from  the 
perfect  conductor  case.  As  such,  the  power  density  calcu¬ 
lated  represents  only  a  first  approximation.  Numerical 
c a k ulat ions ,  which  are  believed  to  be  the  first  made  for 
the  geometry  and  conditions  mentioned  above,  are  presented. 

In  addition  to  the  wedge,  two  other  geometries  considered 
are  the  cone  (representing  waxicons)  and  the  parabolic 
cylinder  'representing  a  flat  surface  with  a  rounded  edge). 
Numerical  computations  using  the  computer  are  difficult  to 
obtain  for  these  two  geometries,  therefore,  only  the 
background  theory  is  presented. 


This  thesis  consists  of  six  chapters  and  four 
appendices.  The  chapter  immediately  following  this 
introduction  presents  the  basic  theory  necessary  to  solve 
the  wave  equation  for  the  three  geometries  specified  in 
the  following  chapter.  Chapter  III.  Chapter  IV  establishes 
criteria  to  verify  computer  routines  to  obtain  the 
numerical  results  for  the  wedge  presented  in  Chapter  V. 

The  final  chapter  presents  the  conclusions  that  can  be 
drawn  from  Chapter  V  and  recommendations  for  further 
study  in  this  area.  The  four  appendices  are  in  support 
of  the  thesis. 
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I I .  General  Background  Theory 

The  purpose  of  this  chapter  is  to  present  the  general 
theory  necessary  to  solve  the  wave  equations  for  the  specific 
geometries  discussed  in  the  next  chapter.  This  is  accomplished 
by  first  presenting  Maxwell's  equations  describing  the  electric 
and  magnetic  fields  (I?  and  ft)  for  all  space  and  time.  The 
geometry  is  that  of  a  perfect  electric  conductor  imbedded  in 
free  space  and  excited  by  a  monochromatic,  time-harmonic 
electromagnetic  wave.  Once  E  and  H  have  been  found,  the 
boundary  conditions  are  used  to  find  the  induced  surface 
current  density  (Jg}  .  Finally,  an  approximation  to 
determine  the  time-averaged  power  density  dissipated  on  the 
surface  of  a  good  conductor  is  derived. 

Maxwell's  Equations  and  Boundary  Conditions 

4 

In  reference  to  Figure  1,  the  following  are  the  assump¬ 
tions  to  be  made:  (1)  Region  I  is  assumed  to  be  free  space 

with  a  conductivity  cr  =  o  mhos/meter  ,  a  permittivity 
- 1 2 

£c  =  8.854  x  10  farad/meter,  and  a  permeability 

__  7 

/^o=  4TT  x  10  henry/meter;  (2)  Region  II  is  assumed  to 
be  a  perfect  electric  conductor  with  conductivity  <T  =  oo  t 
a  permittivity  ,  and  a  permeability  £  ;  (3)  the 

source  is  assumed  to  be  an  infinite  distance  away  from 
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Region  II  and  emitting  a  monochromatic,  time-harmonic  electro¬ 
magnetic  (E-M)  wave.  The  E-M  wave  is  assumed  to  have  elWt 
time  dependence,  where  CU  is  the  radian  frequency  (in 
radians/second)  of  the  E-M  wave.  If  the  incident  electric 
field  E  (where  the  superscript  i  represents  incident) 
is  polarized  such  that  no  component  of  E  is  tangential 
to  the  conductor  surface,  then  the  incident  E-M  is  said  to 

be  transverse  magnetic  (TM) .  If  the  incident  magnetic  field 

^  i 

i!  is  polarized  such  that  no  component  of  H  is  tangential 
to  the  conductor  surface,  then  the  incident  E-M  wave  is  said 
to  be  transverse  electric  (TE) .  The  TM  case  is  considered 
in  the  derivations  to  follow,  unless  stated  otherwise. 


.f  1.  Perfect  Conductor  Illuminate  I  ir,  Jouree 


Maxwell's  equations  for  free  space  near  the  vicinity 
cf  the  conductor  (hence,  neglecting  the  infinitely  distant 
source)  may  be  written  in  terms  of  E  and  H  as  (Ref  1:37,38) 

V  X  H  =  CD 

V  *  E  =  it  (2> 

V  *  E  =  0  <j> 

V  •  H  =  0  C4, 

Since  E  and  H  have  eiwt  time  dependence  and 

Jk 

assuming  that  E  and  H  can  each  be  expressed  as  the 
product  of  two  functions  such  that  one  function  is  a 
function  of  space  only  and  the  other  is  a  function  of  time 
only,  then  =  iu>E  and  =  iu>II  .  By  also  using  the 
relation  for  free  space  only  (Ref  1:38) 


k  =  C5) 

where  k  is  known  as  the  propagation  constant, 
is  the  inverse  of  the  velocity  of  E-M  waves  in  free  space, 
and  X  is  the  wavelength  of  the  E-M  wave  in  free  space,  then 
Eqs  (1)  and  (2)  may  be  rewritten  as 
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V  *  H  =  i  kg  E 


VxE  s  -ikj%  H 


By  taking  the  curl  o£  Eq  (2')  and  substituting  Eq  (1') 

--  A  ^ 

for  V  ^  H  in  terms  of  E  ,  then 

v*v*£  -  kxE =  o 

is  obtained.  Using  the  vecotr  identity  (Ref  1:38) 

VxV*A  «  V(V*A)-VZA 


and  recalling  that  y*  E  =  0  from  Eq  (3),  then  Eq  (6) 
becomes 


vlE  +  kzE  -  0 

__2. 

Vvhere  y  is  the  Laplacian  operator.  In  a  similar  fashion, 
starting  with  Eq  (1')  leads  to  an  equation  the  same  as 

-A 

Eq  (7),  except  with  E  replaced  by  H  .  Equation  (7) 
is  the  complex  vector  wave  equation  that  must  be  solved  in 
order  to  determine  E  .  Once  E  is  found,  the  substitution 

A  «4  afc 

of  E  into  Eq  (2')  determines  H  .  This  gives  E  and  H 
for  Region  I.  Since  Region  II  is  a  perfect  conductor, 
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J  _» 

E  -  H  =  J  -•  0  in  Region  II,  where  J  represents  the 
current  density  (Ref  3:17). 

To  solve  Eq  (7)  for  all  space,  the  boundary  conditions 
for  the  surface  separating  Region  I  from  Region  II  must  be 
given.  The  boundary  conditions,  referring  to  Figure  1, 
which  must  be  imposed  on  the  surface  of  the  conductor  are 
(Ref  1:34) 


n  x  H 


=  J 


(8) 


A  ? 

n  x  E 


(9) 


where  n  denotes  the  unit  vector  normal  to  the  conductor 
surface  and  pointing  into  Region  I,  and  ?s  denotes  the 
surface  current  density  induced  on  the  conductor  surface 
by  the  H  field. 

Examination  of  Eq  (9)  shows  that  for  any  n  , 

A  ^  J* 

n  x  E  =  (J  .  Thus,  E  must  be  equal  to  zero  except 

A  /  > 

when  E  is  normal  to  the  surface,  in  which  case  n  x  E  would 
necessarily  be  zero.  Thus,  there  is  no  tangential  component 

■Jk  A 

of  E  on  the  conductor  surface.  For  Jg  not  equal  to 
zero,  Eq  (8)  implies  that  components  of  H  normal  to  the 
conductor  surface  are  zero,  but  the  tangential  component 
of  H  on  the  conductor  surface  induces  a  surface  current 

A  ^ 

density  normal  to  both  n  and  H 
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r 

a  jk 

Recalling  that  E  and  H  are  separable  functions  of 
space  and  time,  Eq  (7)  can  be  expressed  as  a  differential 
equation  of  only  spatial  coordinates,  since  the  time 

-X  _» 

dependence  will  cancel  out.  For  example,  if  E  ■  Ej^  (u,v,w)f  (t) 

A 

where  E^UjV.w)  is  a  vector  dependent  on  spatial  coordinates 
u,  v,  and  w,  and  f(t)  is  a  scalar  function  of  time  t  only, 
then  Eq  (7)  becomes 


f (t)  V Ej  (u,v,w)  +  f  (t)k2E1  (u,v,w)  =  0 

For  f(t)  f  0  ,  the  above  equation  becomes 

V^Cu.v.w)  +  k2E1(u,v,w)  =  0  (10) 

Equation  (10)  has  the  form  of  a  homogeneous,  linear 
partial  differential  equation  of  order  two.  It  is  also 
assumed  that  u,‘v,  and  w  are  independent  coordinates  of  a 
general  three-dimensional  orthogonal  coordinate  system  and 
L^fu.v.w)  is  such  that 

Ej  (u  ,v  ,w)  =  Ea(u)Eb(v)Ec(w)  (11) 


whe 


re  E^fu.VjW)  is  the  scalar  function  of  E1(u,v,w)  and 


Ti  (u) 
a  v 


Eb(v)  ,  and  Ec (w)  are  separate  scalar  functions 
of  only  one  variable.  Substitution  of  Eq  (11)  into  Eq  (10) 
then  results  in  being  able  to  express  one  side  of  the  equation 
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depending  on  only  one  variable,  and  the  other  side  of  the 
equation  depending  on  the  other  two  variables.  Thus,  both 
sides  of  the  equation  must  be  equal  to  a  constant  for  any 
u  ,.  v  ,  and  w  .  The  process  is  repeated  until  three 
separate  ordinary  differential  equations  are  obtained  for 
E  (u)  ,  E,  (v)  ,  and  E  (w)  .  Solving  the  differential 

equations  for  E  (u)  ,  E  (v)  ,  and  E  (w)  and  super- 

position  of  the  solutions  in  Eq  (11)  into  Eq  (10)  yields 
the  general  solution  to  E^(u,v,w) 

The  method  just  described  is  called  the  method  of 
Separation  of  Variables  (Ref  2:260).  The  coordinate 
systems  for  the  specific  geometries  discussed  in  the  next 
chapter  are  such  that  Eq  (7)  can  be  solved  by  the  method 
of  Separation  of  Variables. 

Finite  Conductivity  Approximations  and  Scattering 

As  stated  ill  the  Introduction,  the  main  objective  of 
this  thesis  is  to  find  an  approixmat ion  to  the  time-averaged 
power  density  dissipated  on  the  surface  of  a  wedge  of  high 
conductivity.  This  is  accomplished  by  using  a  parameter 
known  as  the  skin  depth  and  the  surface  current  density 
produced  on  the  surface  of  a  perfect  conductor  as  in  Eq  (8). 
Because  of  these  two  approximations,  the  following  deriva¬ 
tions  presented  for  the  time -averaged  power  density  loss  is 
just  a  first  approximation  to  the  actual  value. 
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In  reference  to  Figure  1,  it  is  now  assumed  that  the 
conductor  of  Region  II  has  a  finite  conductivity,  O"  . 

It  is  further  assumed  that  Region  II  is  a  good  conductor, 
i.e.,  cr»0O€  .  As  stated  in  the  Introduction,  the 

incident  radiation  of  interest  is  in  the  micron  wavelength 

-  6  14 

(10  meters)  range,  hence  tO  will  be  on  the  order  of  10 

radians  per  second.  The  permittivity  £  of  a  good  conductor 

is  hard  to  measure,  but  is  generally  taken  to  be  on  the  order 

of  £0  for  free  space  (Ref  1:6).  Thus,  for  copper  of 

J 

CT  =  5.8  x  10  mhos  per  meter,  the  condition  ff  »6t)£  is 
well  satisfied.  It  is  also  assumed  that  the  conductor  is 
a  linear  isotropic  medium. 

For  good  conductors,  the  surface  field  amplitude  is 
attenuated  by  a  factor  of  e*1  (  0.368)  after  it  has 

penetrated  a  distance 


into  *  he  conductor  (Ref  1:53).  The  symbol  3  is  known  as 
the  skin  depth  and  is  the  wavelength  of  the  E-M  wave 

in  the  conductor. 

In  reference  to  Figure  2,  if  a  plane  wave  is  incident 
on  the  surface  normal ,  then  Snell's  Law  may  be  written  in 
the  form  (Ref  1:58) 
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The  wave  number,  k  ,  is  a  complex  number  in  general 
It  k  is  then  defined  to  be  (Ref  1:48) 

k  =  k^  -  i  kj 

where  ,  called  the  intrinsic  phase  constant,  is  the 

re.il  part  of  k  ,  and  k^  ,  called  the  intrinsic  attenua 
t i  m  constant,  is  the  imaginary  part  of  k  ,  then  for  a 


good  conductor,  k^  =  =  ( U)JA  CT/  2  )'L  (Ref  1:50).  From 

Eq  (5)  it  is  seen  that  for  free  space,  kR  =  u) 
and  kj.  =  0  .  Eq  (13)  is  now  written  as 


sin  9t  ,  [lQ£0/J0'  f  \ 

sin  0j  "y  O/X  ) 

Since  the  permeability  /*  of  a  linear  medium  is 
approximately  equal  to  /*•  of  free  space  (Ref  1:6)  and 
noting  (i  +  i)/^T*  ,  then  Snell’s  Law  for  a  free 

space  -  good  conductor  interface  is  approximated  as 


n 

Since  for  a  good  conductor  CT»LO£  (for  <T  of  order  10  , 
CO  of  order  10*** ,  and  £0  of  order  10  *^),  then  Eq  (16) 
implies  that  the  ratio  of  sin  0^  to  sin  0-  will  be 
very  small.  For  @.  limited  between  0  and  1T7  2  radians, 
Qj.  must  approximately  equal  zero  for  arbitrary  0^ 

Thus,  the  rapidly  attenuated  surface  fields  propagate  very 
nearly  normal  (with  respect  to  the  surface  tangent)  into  the 
conductor . 

The  surface  current  density  in  Eq  (8)  thus  produces 
a  tangential  electric  field  given  by  (Ref  3:18) 


A 

n  x 


"tan 


c 


> 

J 


(17) 


13 


It  is  useful  to  use  the  concept  of  scattering  in  order 
to  find  H  everywhere  so  Eq  (2')  can  he  used  to  find  H 
for  use  in  Eq  (18).  The  scattered  field  can  be  thought  of 
as  the  field  produced  by  the  currents  induced  on  the 
conductor  surface  as  in  Eq  (8).  This  scattered  field  (Es  ) 
is  thus  equal  to  the  difference  between  the  field  everywhere 
with  the  conductor  present  (T:)  and  the  incident  field  with 
the  conductor  absent  (E  )  .  This  is  written  as  (Refl:113) 


14 


(19) 


Jkc  A  ][  «Jb 

Thus,  by  solving  for  E  and  given  E  ,  E  is  found  and 
Eq  (2')  then  yields  U  .  The  solution  for  H  is  then  used 
in  Eq  (18)  to  give  a  first  approximation  of  the  time -averaged 
power  density  dissipated  on  the  highly  conductive  conductor 
surface . 

The  necessary  background  theory  has  now  been  presented 
to  successfully  derive  the  solutions  to  the  specific  geometries 
presented  in  the  next  chapter.  Since  many  of  the  sharp  edges 
occur  in  highly  conductive  plane  (flat  surface)  mirrors  (see 
Figure  3),  the  first  geometry  discussed  in  the  next  chapter 
is  that  of  an  infinite  wedge,  which  has  two  flat  surfaces 
meeting  at  a  sharp  edge.  Since  present  technology  is  such 
that  perfectly  sharp  edges  are  not  possible,  but  in  fact  the 
edges  are  rounded,  then  the  geometry  of  a  parabolic -cyl inder 
is  used  to  approximate  flat  conductor  surfaces  with  rounded 
corners.  Finally,  since  waxicons  (see  Figure  4)  are  being 
preposed  for  use  in  the  resonators  of  high  power  chemical 
lasers,  the  cone  is  used  to  approximate  such  optical 
components . 
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III.  Theory  Extended  to  Specific  Geometries 

•  The  purpose  of  this  chapter  is  to  specialize  the 
theory  in  the  last  chapter  for  specific  boundary  conditions. 
The  boundary  conditions  lead  to  the  choice  of  what  coordinate 
system  is  used.  The  first  section  considers  the  case  of 
scattering  by  an  infinite,  perfectly  conducting  wedge 
illuminated  by  a  parallel  electric  line  source.  The  line 
source  is  then  moved  to  infinity  so  as  to  approximate  plane 
wave  incidence.  A  first  approximation  is  then  made  for 
determining  the  time-averaged  power  density  dissipated  on 
the  surface  of  a  highly  conductive  wedge.  As  a  special  case 
of  the  wedge,  the  half -plane  problem  using  the  wedge  solution 
is  compared  to  the  classical  solution  for  the  half-plane 
(horn  and  Wolf’s  derivation  is  used;  Ref  6:559-570).  The 
method  of  expanding  the  infinite  wedge  problem  to  one  of  a 
finite  wedge  is  introduced. 

The  second  section  considers  the  wave  equation  for  a 
parabolic -cylindrical  coordinate  system.  Representing  a 
bend  with  a  rounded  corner,  the  parabolic  cylinder  is 
discussed;  in  particular,  the  difficulty  in  obtaining 
numerical  results  to  the  solution  of  the  parabol ic -cylindrical 
wave  equation  is  identified.  Finally,  the  last  section 
considers  the  wave  equation  for  spherical  coordinates.  A 


17 


f 


brief  discussion  to  the  solution  of  the  spherical  wave 
equation  for  a  cone  boundary  (cones  have  boundaries  easily 
expressed  in  spherical  coordinates)  is  discussed. 

The  first  section  is  the  most  thoroughly  discussed 
section,  since  the  solutions  derived  are  the  basis  for 
the  next  chapter  and  for  the  numerical  results  presented 
ir>  Chapter  V. 

C  1  indrical  Coordinates  and  the  hedge 


The  wedge  has  surfaces  that  are  conveniently  expressed 
in  cylindrical  coordinates  (see  Figure  5b).  The  cylindrical 
u.i  o  equation  is  introduced  first,  followed  by  scattering 
by  wedges.  The  cylindrical  coordinate  system  used  is 
ck  l  ined  in  Figure  5a. 


Cylindrical  Wave  Equation.  The  scalar  version  of 
Eq  (7)  in  cylindrical  coordinates  with  the  Laplacian 
e\j  vessed  in  cylindrical  coordinates  and  £  replaced  by 
t!..  scalar  function  ^  is  (Ref  1:198) 


f  r  f 


(20) 


Eq  (20)  is  also  known  as  the  scalar  Helmholtz  equation  in 
cylindrical  coordinates.  If  is  separable  into  the  form 
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Figure  5b.  diedge  as  Special  Case  of  Cylinder 


'/'((.ti  z)=  'f=  R(f)§(^)Z(z)*R^ 

and  1 4  (20)  is  assumed  to  have  e1U>t  time  dependence,  the 
method  of  Separation  of  Variables  can  be  used.  The  well 
known  separable  equations  resulting  from  substitution  of 
Eq  (21)  into  Eq  (20)  are 


/_§ 

u 


j  +  n2  £  ~  0 


fi(e#MM'-"2]R’0 


.chore  k  in  Eq  (24)  is  defined  by 


e 


k2  +  k2 

e 


(23) 

(24) 


(25) 


Equations  (22)  and  (23)  are  homogeneous,  linear 
oie inary  differential  equations  with  constant  coefficients, 
iii.-o  give  rise  to  harmonic  solutions  (sines  and  cosines), 
iui.ation  (24)  is  Bessel's  equation  of  order  n  and  has 
Be.-sel  function  solutions  (Ref  1:199).  Equation  (21)  then 
becomes 


z)  -  Bn  (kff)h(n^)h(kzz) 


(26) 


where  h(n^)  and  h(kzz)  represent  harmonic  functions  for 
tirv  solutions  of  Eqs  (23)  and  (22),  respectively.  Bn(kf^>) 
represents  a  Bessel  function  solution  to  Eq  (24)  (see 
Appendix  A  for  the  definitions  of  various  Bessel  functions). 
From  the  principle  of  superposition,  linear  combinations  of 
Eq  (26)  are  also  solutions  to  Eq  (20).  The  possible  values 
for  n  ,  k  ,  and  k  are  called  eigenvalues  and  any 

r  2 

two  may  be  summed  or  integrated  over,  except  k  and  k 

f 
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since  they  are  interrelated  by  Eq  (25).  Equation  (26)  is 
specified  in  the  discussion  to  follow. 

Scattering  by  a  h'edge.  Consider  the  case  of  an  electric 
current  filament  of  current  1  (in  amps)  located  parallel 
to  a  perfectly  conducting  wedge,  as  shown  in  Figure  6. 

The  wedge  has  an  angle  2e(  and  covers  the  space  (in 
cylindrical  coordinates  ,  <f>  ,  z  ):  2TT'-  oi  ^  (J>  ^  & (  > 

0  ^  p  <  oo  ,  and  —  oo< z<  oo  •  The  filament  and  wedge  are 
assumed  to  be  imbedded  in  free  space  (CT=  0,  e  '£„v“»A>  ■ 
The  wedge  is  assumed  to  be  a  linear,  isotropic  perfect 
conductor  (CT=©o)  .  The  incident  electric  field  can  be 

described  as  an  outward-traveling  cylindrical  wave,  hence 
the  Hankel  function  of  the  second  kind  and  order  n  (n  an 
integer),  denoted  by  H^2^(k  )  ,  is  used  (provided  eiCtft 

t ime -dependence  is  assumed)  (Ref  1:201).  The  incident  electric 


6.  perfectly  conducting  v 1U;  u-  urrer,t  filament 


:*i 


1 


A 


(27) 


i  ield  is  z-polarized  and  given  by  (Ref  1:236) 

Ei  -  HfaiH'l) 

v.r.ere  -  £.  ^  ^  ^  “  2.  ^  ^  C.OS  and  is 

evident  from  the  geometry  in  Figure  6. 

The  addition  theorem  for  Hankel  functions  is  given 

7-  e<e' 

.e  Jn( and  Jn(k^)  are  Bessel  functions  of  the 
i  >t  kind  of  integer  order  n  .  Equation  (28)  substituted 
ir.’.o  Eq  (27)  then  gives  the  incident  field  as 


e:-< 


-S  f>f. 


(29) 


The  total  electric  field  everywhere,  Ez  ,  is  given 
by  Eq  (19)  as 


E  =  E1  ♦  ES 
z  z  z 


(30) 
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is  considered  as 


r.ere  E*  is  given  by  Hq  (2  91  ann  ir 

t  ,'c  field  originating  from  the  induced  currents  on  the 

surface  of  the  wedge  and  is  called  the  scattered  field 

(Kef  t>:SS8).  The  scattered  electric  field  is  in  the  form 

of  hq  (29),  but  with  the  integer  order  n  changed  to  real 

(21 

oiler  v  and  replaced  by  II  (k^)  .  The  order 

v  is  determined  (to  follow  laterl  by  the  boundary  conditions 

on  the  wedge  surfaces.  The  scattered  electric  field  E^ 

then  added  to  the  incident  electric  field  E*  to  yield 

tie  total  field  ,  as  in  Eq  (30).  Since  the  boundary 

condition  must  have  =  E1  at  the  wedge  surfaces,  the 

z  z 

i  i  ..  I  field  must  also  be  z -polarized  (Ref  1:238  ).  The 
total  field  is  thus  written  as  (Ref  1:238) 


I  QvHr(ke')Xfrf)sinLv(f-«o>in&w^ 
,(2) 


Ez=  \ 


(31) 


z  a  J>e')  TO  rf' 


1  he  boundary  conditions  t  ha  t  Miust  he  imposed  on  the 
faces  of  the  perfectly  conducting  wedge  are 


Ez  -  0  ,  ^=c* 

ez  =  o  ,  2ir-o< 


(32) 
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.‘iolving  Hq  (31)  for  the  above  boundary  conditions  implies 


sin 

[v(^-«0]  =  0  ,  $ 

=  * 

(33a) 

s  in 

'’Q. 

at 

o 

II 

/—X 

« 

1 

> 

=  air-  * 

(33b) 

quat ion 

(33a)  is  true  for  any 

v  ,  but  Eq  (33b) 

is  true 

on i y  f or 


v(21T-0<-00  =  2v(fT-«n  =  mlT  (34) 

eh.  re  m  is  an  integer.  The  value  m  =  0  (thus  v  =  0  ) 
is  not  allowed  because  E  in  hq  (31)  becomes  zero  for  any 
(  ^)  .  The  values  of  m  equal  to  a  negative  integer  are 

jImi  not  allowed,  since  by  the  definition  of  Jy(kp)  for 
V  negative  (see  Appendix  A),  li  in  liq  (31)  tends  to 
ii.iiuity  as  k  tends  to  zero.  Thus,  liq  (34)  gives  the 
v  .  .  ;i v  ,t  1  ues  v  of  liq  (31)  as 

mqTT 

v  ‘  2 (ir- a)  ’  m  3'  1351 

The  constants  ay  are  determined  by  the  nature  of 
the  source  (Ref  1:239).  In  reference  to  Figure  7a,  the  current 
filament  is  considered  to  be  an  impulse  of  current  of  strength 
I  (in  amps).  Figure  7b  represents  a  vanishingly  small 


,ii  . 'a  surrounding  the  current  filam-nt;  small  enough  so  that 
ii,e  current  filament  can  be  considered  as  a  current  sheet 
i  ^  *  p* )  separating  the  regions  for  p  >  p'  and  £  <  q*  . 
\ppendix  B  provides  the  details  which,  by  referring  to 
1  igure  7b,  leads  to  the  source  current  density  as 


J 


z 


— - ;  Va  sinjvl  ^'“CX)  ]sin[v(^-C<  )]  (36) 

7  v 


equation  (36)  is  in  the  form  of  a  Fourier  sine  series  for 
the  current  on  ^  ^  .  The  Fourier  sine  series  for  an 
impulse  of  current  of  strength  I  (in  amps)  at  e  -  e' 

.  i  i  i  j  i'  is  given  by  (Ref  1:23d) 


.1 

z 


X 

(ir-cty 


sin  [v(jT-  cx  )  ]  sin  fvfjtf-  <X  )  ] 

V 


(37) 


'..■  mpai  ing  Fq  (36)  with  Iiq  (37)  thus  yields 

OJ//.TT  I 
2(tr-o<) 


(38) 


Of  special  interest  to  this  tliesis  is  the  case  of 

tin  nochromatic  plane-wave  illumination.  This  is  accomplished 

!m  letting  the  current  filament  recede  to  infinity,  i.e., 

kp"* _ >oO  .  The  incident  electric  field  in  Fq  (27)  thus 

( 2 ) 

in  omes  (using  the  large  argument  approximation  for  lr  (x) 

,i  x — >oo  given  in  Appendix  A)  (Ref  1:240) 


c6 


(39) 


rl 


^HTUU-rl) 

O yu0I 


€L  6, 


I  quation  (39)  can  thus  be  written  in  the  form  of  a  plane 
iiave  as  (Ref  1:240) 


l 


z 


E  eil<ecos(*-^ 

o 


(40) 


Only  the  equation  for  £  <  need  be  considered  in 

!.q  (31)  to  find  the  electric  field  !•„  everywhere  for 

k  - >QO‘  Using  the  asymptotic  approximation  for  H^^(k^') 

given  in  Appendix  A  as  (Ref  1:463) 


( A  - 1 2 ) 


>1  (31)  then  becomes 


Substitution  of  ay  given  by  l.q  (38)  and  nQ  given  by 
q  (4  0)  into  F.q  (41)  yields 


2? 


Ez  =  2iV  Jv(kP1s5n!v(^"°<)lsinlv(^"0(  )J  (42) 

Tr- oC  y 


where  v  is  given  by  Eq  (35).  Equation  (42)  thus  represents 
the  solution  of  the  electric  field  for  a  monochromatic, 
t-polarized  plane  wave  incident  on  a  perfectly  conducting 
wedge  of  angle  2  (see  Figure  6) .  The  solution  for  the 
magnetic  field  H  is  obtained  by  applying  Maxwell's  equation, 
I q  (2  '  )  ,  to  Eq  (42) ,  i. e. , 


S  .  i  /  *• '  (  A  3  Ex  A  _  'J  ) 

T?a(  50  *) 


! ho  surface  current  density,  J  ,  is  obtained  by  substitu¬ 
ting  Eq  (43)  into 


A  -fc 

n  x  H 


or  the  top  surface  of  the  wedge,  n  =  a- 


The  a. 


«  v  j  t  t  iiu  i  i  i  uv.  v  i  i  m;  (  u  ~  a  ^  •  i  lit'  <s  ^ 

component  of  H  in  Eq  (43)  is  the  only  term  used  since 
Vj  K  a^  =  0  .  The  surface  magnetic  field,  denoted  by  Hs 

is  thus  given  as 


A  [JZ  i_  n 

V.-J  M*  P  id  f 


where,  after  differentiating  hq  (42)  with  respect  to  $  , 
the  boundary  condition  was  imposed.  The  surface 

current  density,  Jg  ,  is  thus  given  by 


-»  /A  A 

J  =  J  a  =  -  H  a„ 

S  Z  Z  fs  2 


u  here 


If  the  wedge  is  a  good  conductor  of  finite  conductivity 
such  that  CT  »co£  ,  then  hq  (18)  of  Chapter  II  is  used 
to  obtain  a  first  approximation  to  the  time-averaged  power 
density  dissipated  on  the  surface  of  the  wedge,  provided 
that  the  Leontovich  Boundary  Condition  is  satisfied  (see 
Appendix  B) .  It  is  assumed  throughout  the  rest  of  this 
thesis  that  the  Leontovich  Boundary  Condition  is  satisfied. 
Recalling  Lqs  (12)  and  (18) 


*“  (■ 


2  \y2 


°  ~  \u)jua  j 

^.„= 

where  |lIo|  =  t  lie  dot  product  of 

and  its  complex  conjugate  given  by  l:.q  (45),  thus  gives 
the  time-averaged  power  density  loss  on  the  top  wedge 
surface  as 
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(46) 


F  =  4-f^Y- 

Ui«“  <*d/io  Up  l*-*)  1 


£  i v V  X  f)  sin  [v^-o<j] 

V 

2(-i)VvX(kf)sin[v(^][ 


The  problem  for  the  wedge  just  discussed  is  for  the 
I'M  case.  The  TE  case  is  solved  in  a  similar  manner  as  the 
1M  case,  using  a  magnetic  current  filament  in  place  of  the 
electric  one.  For  plane  wave  incidence,  the  result  for  the 
magnetic  field  is  given  by  (Ref  1:242) 


H 


ttM 

v(k^)cos  [v(jT-0<  )]cos[v(j£-c()] 
1,  v  *  0 


where 


£v  = 


and 


2,  v  >  0 


mir 


(47) 


v 


2  (IT- o£)  * 


III  =  0,  1,  2, 


At  the  top  wedge  surface  )  ,  t  iie  resulting  magnetic 

field  on  the  surface,  ii  ,  becomes 


IT  Ho 


"ZS  =  22 €viVjv ck^) cos  l V (^—  04 )  ] 


(48) 


Substitution  of  Eq  (48)  into  Eq  (10)  yields  the  induced 


surface  current  density  Jg  as 


ir  M 

js  =  Hz  ■  “fiizr  Eeviv.Jvlk^)c°stv(^/-00]  (49) 

liquations  (12)  and  (18)  are  used  with  Rq  (48)  to  obtain 
a  first  approximation  for  the  time-averaged  power  density 
dissipated  on  the  top  surface  of  the  highly  conductive 
wedge.  Of  note  is  the  differences  in  the  TM  case  of  Eq  (44) 
to  the  TE  case  of  Eq  (48).  The  TM  case  of  Eq  (44)  has  the 
term  kp  in  the  denominator  while  the  TE  case  of  Eq  (48) 
does  not.  Equation  (44)  also  varies  with  sin  [v(^/-oC  )  ] 
while  Eq  (48)  varies  with  cos  [vf^'-c*  )  ]  .  In  addition, 

!  q  148)  has  an  additional  v  =  0  term  and  slightly  differ¬ 
ent  constants. 

The  TE  case  is  mentioned  here  only  to  present  the 
equations  necessary  to  determine  the  time-averaged  power 
density  loss  due  to  TE  illumination  so  as  to  compare  to 
the  TM  results  presented  in  Chapter  V. 

Hal  f-Pl ane .  For  o(  =  0  ,  t  lie  wedge  in  Figure  6 

reduces  to  the  classical  half-plane  problem.  The  half¬ 
plane  perfect  conductor  covers  the  space:  0  ^  x  <  ©o  , 
y  =  0  ,  and  -oo<z<00  ;  see  Figure  8. 

It  is  assumed,  in  reference  to  Figure  8,  that  a 
monochromatic  plane  wave  (the  electric  field  ^-polarized) 
is  incident  at  an  angle  on  a  perfectly  conducting 

half-plane  imbedded  in  free  space.  The  incident  electric 
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figure  8.  Perfectly  conducting  half-plane  illuminated 
by  z-polarized  E  and  -polarized  H. 


field  is  given  by  Eq  (40).  The  total  electric  field  every¬ 
where  is  then  given  by  Eq  (42)  with  o(  =  0  .  The  eigen¬ 

values  v  given  in  Eq  (35)  reduce  to 


v 


m'fr 

aCtr-oO 


n 


3, 


(50) 


1  he  total  electric  field  is  then  written  as 


2  E 


CO 

„r 


■ 


(51) 


Specializing  Eq  (45)  for  oC  =  0  gives  the  induced 
surface  current  density  as 


3- 


(52) 


An  alternative  solution  for  the  half-plane  is  treated 
thoroughly  by  Born  and  Wolf  (Ref  6:559-584).  The  method 
involves  using  dual  integral  equations  for  the  E  and  H 
fields,  of  which  dual  is  meant  that  a  single  unknown  function 
satisfies  different  equations  for  two  distinct  ranges  of 
parameter  (Ref  6:564-565).  The  final  result  is  Sommerfeld's 
equation  (Ref  6:569) 


-here  and  <*) 


are  the  1  resnel  cosine  and 


sme  integrals,  respectively,  and  defined  as  (Ref  5:300) 
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(54) 


CM-  fcos(ft2)jt 

*'o 

JW-  /"»(! t')jt 

'o 

Equation  (53)  is  valid  for  the  incident  field  E*  given 
hv  Eq  (40).  (Born  and  Wolf  used  =  exp[-ik^  cos  ( -  f>)  \  ) 
and  for  Eq  normalized  to  equal  one. 

By  using  Maxwell's  equation,  Eq  (2')  again,  and  the 
boundary  condition  Eq  (8)  to  get  the  surface  current  density, 
t  hen 


In  the  next  chapter,  a  truncated  version  of  Eq  (52) 
is  compared  numerically  on  the  computer  against  Eq  (55)  so 
as  to  give  a  clue  on  the  reliability  of  the  truncation  and 
Bessel  function  subroutine.  To  get  an  indication  of  the 
reliability  of  the  numerical  results  presented  in  Chapter  V 
and  based  on  Eq  (46)  for  the  wedge,  a  comparison  of 


3d 


For  (X  =  90°,  Fq  (35)  for  v  reduces  to 


Wit 


2(lT-oO 

and  Fq  (45)  for  o(=ir/2  reduces  to 


=  m  ,  m  =  (1,2,3,...)  (56) 


'<*! 


oo 


Using  cartesian  coordinates,  the  classical  solution  to 


tins  problem  is  (Ref  16:473,  531) 


C  0  S  ^  0; 


(S8) 


where  9«  is  the  angle  of  incidence  with  respect  to  the 
surface  normal  (  e,  -  4'  -  fr )  .  Comparison  of  Fq  (58) 


with  Fq  (57)  implies  H  -II,  (evident  from  Figure  9  is 

'  s  ^  s 

that  a,  is  in  the  a  direction  on  the  conductor  surface) 

y 


A  truncated  version  of  Fq  (57)  is  compared  numerically  to 
lq  (58)  in  the  next  chapter. 


extension  to  Three  Dimensions.  Now  consider  the  wedge 
in  Figure  10  that  covers  the  space:  0^  ^  <  oo  f 

z  ^  ^  z  z7  .  It  has  the  same  g comet ry 
as  the  wedge  in  Figure  6,  except  that  now  the  wedge  is 
finite  along  the  z-axis  .  Thus,  a  three-dimensional  wave 


equation  must  be  solved  instead  of  the  simpler  two-dimensional 
wave  equation  for  the  case  of  non-varying  geometry  with 
respect  to  the  z-axis.  However,  by  applying  a  Fourier 
transformation  with  respect  to  the  z-axis,  the  three-dimen¬ 
sional  wave  equation  can  be  reduced  to  a  two-dimensional  wave 
equation. 


As  an  example  of  this, 
the  three-dimensional  wave 


if  f ({,#,-) 

equat  ion 


is  a  solution  to 


tt  *  k  1t-  o 


then  ¥<(.*»>  ■  fti t.f,  z )  e  '  h  “  d  z 

—  OO 

the  two-dimensional  wave  equation 


is  a  solution  to 
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V*  ¥  +  K*  f  -  o 


K2  2  2 

•  k  -  w  (Ref  1:243).  When  the  two-dimensional 
have  equation  above  is  solved  for  ,  then  f'  may  be 

solved  using  the  inversion  (Ref  1:243) 


,  iwz  . 
w)  e  dw 


(59) 


-  o© 


Parabolic  Cylinder 

Since  technology  has  not  advanced  yet  to  the  stage  of 
producing  perfectly  sharp  edges,  the  parabolic  cylinder  is 
useful  to  approximate  a  rounded  corner.  The  two-dimensional 
parabolic-cylindrical  coordinate  system  is  shown  in  Figure  11. 
Thu  z-axis  is  perpendicular  to  the  plane  of  the  paper.  The 
transformation  from  rectangular  coordinates  (x,y,z)  to 
;  .i  rabolic-cyl indrical  coordinates  (u,v,z)  is  given  in  Figure  11 

The  wave  equation  that  must  be  solved,  Fq  (7),  in  parabolic 
cylindrical  coordinates  is  (Ref  7:282) 


L  (t) 

c2(u2fvz)  (  3u 


tt  ,  ?t\  ,  ft .  i2. 


+  — — r  I  +  — ,  +k  t  =  0 


2  z 


(60) 


■a  no  re  ( u ,  v ,  z )  =  f*  replaces  the  vector  "E  .  Assuming 

,  t j 

f iu,v,z)  is  separable,  i.e.. 
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Pi^ure  11.  Parabolic  cylinder  at  v  -  vQ 


f(u,v,z)  =  U(u)V(v)Z(z) 


(61) 


i hen  the  mot  hod  of  Separation  of  Variables  yields  the  three 
differential  equations  (see  Appendix  11  lor  details) 


&  „  £  Z  -o 

cl  z 

+  ( z  v  4- 1  -  u2)  U  =  0 

Ju1  V  ' 


(63) 


i^L  +  [-(av  +  1)-  v*]V  =  0 


(64) 


liquation  (62)  is  the  harmonic  equation  discussed 
earlier,  leading  to  linear  combinations  of  cosines  and 
sines.  These  solutions  are  bounded  Cor  —  o©<  z  <  oo  . 
liquation  (62)  has  as  a  general  solution  (Ref  7:294) 


U(u) 


(65) 


v.iiore  k2  is  given  in  Appendix  B  and  11^(0  is  called 
ihe  liermite  function  of  degree  V  and  defined  in 
Appendix  A. 

The  solution  to  Eq  (65)  must  be  bounded  for  u 
approaching  i  oO  .  Using  the  asymptotic  expansions  for 
ll^(ix)  and  H  (  (ix)  given  in  Appendix  A  thus  forces 
'•  in  Eq  (65)  to  equal  zero  and  V  to  equal  non-negative 
integers.  Thus,  the  condition  of  finite  U(u)  forces 
l.q  (65)  to  become 


U(u) 


-k*c  u  z/Z 

A  e 


HnS/C?U> 


n  =  (0,1 ,2 ,3, ... ) 


(65' 


With  the  restriction  of  V  =  n  =  ( 0 , 1  ,2 ,3  , . . . 


v I  (b-4)  has  the  general  solution  (Ref  7  :  294) 

v i v t .  V£V) 

+Qek*‘vV2Hn.i(Vi?v)  1661 

•  -  ir.g  the  asymptotic  expansions  for  Ilf  i x )  and 
!  given  in  Appendix  A  and  the  fact  that  V ( v ) 

.ust  he  hounded,  yields 


(66') 


Die  surface  of  the  conductor  (see  Figure  11 j  is  at 
v  -  v  •  The  boundary  condition  that  must  be  satisfied 
i  -  t  hen 


Using  superposition  of  the  solutions  to  Eqs  (62), 

( o 3 ) ,  and  (64),  setting  z  =  0  for  the  two-dimensional 
problem,  and  combining  the  constant  A  in  Eq  1 05')  with 
?  lie  constant  Q  in  Eq  (66')  to  make  n 


then  Eq  (67) 


o°  oo  _  _ 

»  ‘  SI-;,, (68» 

ro«o  n*o 

t  is  mM'ssarv  to  find  the  coefficients  0  ,  for  each 

Mn,n 

and  n  that  satisfy  liq  (  o  s )  in  order  to  complete  the 
solution  to  liq  (60).  The  boundary  condition  of  liq  (67) 
forces  the  conditions  of  liq  (('8).  finding  solutions  to 
n  on  a  computer  presents  a  more  difficult  problem 
than  that  of  ho  (45). 

'  iheiical  Coordinates  and  the  Cone 

In  this  section,  the  spherical  wave  equation  is  solved 
,.,d  the  solution  for  the  boundary  conditions  of  a  perfectl 
conducting  cone  are  briefly  discussed.  The  spherical 
coordinate  system  to  be  used,  with  coordinates  (r  , 

is  defined  in  Figure  12. 


Figure  17.  Jpheri.  >1  Coordina  to  Jyctera 


The  scalar  Helmholtz  equation 


Spherica l_  Wave  Equation. 
n  spherical  coordinates  is  given  by  (Ref  1:264) 


0= 


4-  —  (c 

r  ir  V 


4- 


r2Sin0  ^0 


Tr 

L  -ff 

*  f2,  si Y]z&  0  ^ 

is  sum  ing  ft  r,0,jzJ)  is  separable,  i.e. 


(si"9lS 


(69) 


-  f  =  R(r)0(e)i(^)=  R©$ 


( b9'  ) 


then  the  method  of  Separation  of  Variables  yields  the  three 
l  i  f  l'e  i  ential  equations  (Ref  1:265) 


/j 


+  m2  §  =0 


jr  (rZ~^r)  +  [kV2  “  n(n4'13R  =  0 

7e(si"e#)  +  KHl)'ym4]e’° 


(70) 


(71) 


(  "2  ) 


sine  o/e 


liquation  (70)  is  the  harmonic  equation  seen  in  the 
pievious  two  sections.  Its  solutions  are  denoted  by  h(m^)  . 
equation  (71)  is  closely  related  to  Bessel's  equation  as 
in  liq  (24)  of  the  first  section.  It  has  solutions  that 
are  called  spherical  Bessel  functions.  Denoting  the 
in  dinary  Bessel  functions  by  B  (kr)  and  b^fkr)  for 


he  spherical  Bessel  functions,  the  relation  between  the 
v. o  is  given  by  (Ref  1  :  265) 

b"(kr)  *J l^7rwkr) 

quation  (72)  is  related  to  Legendre's  equation  and  gives 
ise  to  associated  Legendre  solutions  (see  Appendix  A), 
.noted  by  I.^(cos0)  (Ref  1  :  265). 

The  solution  to  ^  in  Eq  (69)  is  therefore  written 
s  a  linear  combination  of 

=  h(m^)bn(kr)LjJ(cose)  (73) 

Cone.  Figure  13a  represents  a  finite  perfectly 
inducting  cone  imbedded  in  free  space  and  Figure  13b 
i presents  an  infinite  perfectly  conducting  cone  imbedded 
free  space.  Both  cones  bate  conductor  surfaces  easily 
.,<ressed  in  spherical  coordinates. 

Rogers,  Schindler,  and  Schultz  (Ref  8:67-80)  have 
:  mined  numerical  results  from  a  computer  for  the  case 
f  a  monochromatic  plane  wave  (electric  field  x-polarized) 

,'K  ident  head  on  the  tip  of  the  cone  in  Figure  13a 
,  ropagat ing  in  the  -2  direction).  The  boundary  conditions 
at  must  be  satisfied  include  (Ref  8:70) 


(b)  Infinite  cone  i  1  lnnii.aiv  :  tyj  vii.£  oourte 


I:  =  0  ,  Q  =  0o  and  r  <  b 


uni  finite  energy  at  the  edge  of  the  cone  (r — , 

)-*&)• 

Warrington  (Ref  1:303-300)  examines  the  case  of  a 

A 

in  i  form  jP-directed  ring  source  illuminating  an  infinite 
■  n : : e  as  in  figure  13b.  The  boundary  conditions  that  must 


satisfied  include  (Ref  1:303) 


E 


0 


0=6o 


and  that  the  energy  remain  hounded  ( finite).  The  boundary 
conditions  for  the  cone  in  both  cases  require  that  (Ref  1: 
504  and  Ref  8:72) 

.ole 

.here  I’v(cos  0)  is  called  the  associated  Legendre 
function  of  the  first  kind  of  order  v  (see  Appendix  A 
for  a  definition).  The  important  point  here  is  that  the 
physical  boundaries  of  the  cone  geometry  force  the  condition 
in  Eq  (74)  to  be  true.  Numerical  calculations  are  difficult 
to  obtain  for  E  because  of  the  problem  of  obtaining  the 
eigenvalues  v  and  the  eigenfunctions  Pv(cos0)  in 
1  i  (74)  (Ref  1:305  and  Ref  8:77). 
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IV.  Verification  of  Routines 


Introduction 

The  purpose  of  this  chapter  is  to  verify  a  Bessel 
function  computer  subroutine  able  to  fulfill  the  require¬ 
ments  of  the  criteria  established  in  the  following 
sections.  Equation  (46)  is  the  key  equation  in  support 
of  the  main  objective  of  this  thesis.  ft  involves  fractional 
order  Bessel  functions  of  the  first  kind.  Though  many 
reliable  computer  subroutines  are  available  for  integer 
order  Bessel  functions  and  spherical  (order  of  odd  multiples 
of  one-half)  Bessel  functions  of  the  first  kind,  not  many 
are  available  to  calculate  fractional  order  Bessel  functions. 
Equation  (46)  involves  an  infinite  summation.  Since  com¬ 
puters  can  only  compute  finite  sums,  it  is  necessary  to 
truncate  F.q  (46)  in  such  a  way  that  the  terms  truncated  off 
.no  insignificant  compared  to  the  finite  sum.  Tor  higher 
arguments  of  the  Bessel  function  (on  the  order  of  one 
hundred),  orders  greater  than  the  argument  are  needed  before 
the  value  of  the  Bessel  functions  fall  off  appreciably. 
Therefore,  it  is  also  necessary  to  find  a  computer  subroutine 
that  calculates  high  order  (and  fractional)  Bessel  functions 
of  the  first  kind.  The  availability  of  reliable  computer 
mibrout ines  of  this  kind  is  still  further  reduced.  The 
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necessity  of  finding  such  a  subroutine  is  thus  of  the 
utmost  importance. 

A  thorough  literature  check  has  revealed  no  numerical 
results  for  the  time-averaged  power  density  loss  dissipated 
on  the  surface  of  a  conducting  wedge  near  the  vicinity  of 
its  edge.  Thus,  there  isn’t  any  data  to  check  the  results 
of  a  truncation  of  Eq  (46)  against.  The  next  chapter 
presents  numerical  results  for  a  truncated  version  of 
Eq  (46).  In  order  to  accept  these  results  as  valid, 
alternative  checks  must  be  made.  One  method  involves 
verifying  the  Bessel  function  subroutine  for  the  range  of 
argument  and  order  established  in  the  criteria  of  the  next 
section.  Another  method  involves  the  half-plane  problem 
discussed  in  the  previous  chapter.  Equation  (46),  for 
the  wedge,  and  Eq  (52'),  for  the  half-plane,  are  very 
similar  in  nature.  Both  equations  involve  an  infinite 
sum  containing  a  Bessel  function  of  the  first  kind  of 
increasing  order  (fractional  and/or  integer).  They  are 
also  both  multiplied  by  [  j  j1)  order  #  (order)  .  sin 
l order  .  constant)]  ,  where  the  word  "order"  represents 
the  value  of  the  order.  The  order  is  the  increasing  index 
of  the  summation.  A  finite  summation  version  of  Eq  (52') 
for  the  half-plane  is  compared  against  the  classical 
solution  (from  Born  and  Wolf)  of  Eq  (55').  Thus,  because 
of  the  similarity  in  nature  of  Hqs  (46)  and  (52'),  it  is 
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hoped  that  the  results  of  Eq  f 4 6)  presented  in  the  next 
chapter  can  be  accepted  as  valid. 

Criteria  and  Results  for  Veri  fv ing  Bessel  Function  Subroutine 

This  section  provides  support  to  verify  a  Bessel 

function  subroutine  that  produces  Bessel  functions  of  the 

first  kind,  given  by  J_^^^Ck  ^  )  ,  for  use  in  Eqs  (45) 

and  (46),  where  V  is  the  fractional  part  of  the  order  and 

\  is  the  integer  part  of  the  order. 

The  argument  of  the  Bessel  function  is  ,  where 

K  is  the  wavenumber  given  by  Eq  (5)  in  Chapter  II,  and 

£>  is  the  radial  distance  from  the  edge  to  the  point  of 

interest.  Since  the  main  objective  of  this  thesis  is  to 

find  the  time-averaged  power  density  dissipated  on  the 

surface  of  a  highly  conducting  infinite  wedge  (of  angle 

90°)  near  the  vicinity  of  the  edge,  it  is  necessary  to 

establish  the  range  of  k^  .  Atomic  distances  (radius 

a!  atoms,  bond  lengths  between  atoms  making  up  a  molecule, 

°  -10 

etc.)  range  around  the  order  of  one  angstrom  (1A  =  10 

O 

meters).  The  radius  of  a  copper  ion  (+1)  is  0.96A,  for 
example  (Ref  4:F-171).  Since  t he  incident  radiation  of 
interest  is  a  monochromatic  plane-wave  in  the  micron 
wavelength  range,  a  value  of  k^  =  0,001  represents  a 
point  ^  c*  0.000159  wavelengths  (of  the  incident  radiation) 
from  the  edge.  For  a  CC>2  laser  of  X  =  10.6  microns  , 


r* 


is  Chen  equal  to  approximately  16. 9A  away  from  the 
edge,  which  is  closer  to  the  edge  than  what  present 
technology  can  probably  make  for  a  sharp  corner  of  a  plane, 
flat  mirror.  Increasing  k ^  by  a  factor  of  105  (k^  -  100) 

gives  the  point  ^  —  15.9  wavelengths  away  from  the  edge, 
or  ^  fir  169  microns  =  0.169  mm  for  ^  =  10.6  microns.  Thus, 
by  choosing  the  range  0.001^  kp4  100  ,  an  adequate  range 

is  provided  for  investigating  the  power  density  dissipated 
on  the  surface  of  a  highly  conductive  infinite  wedge  near 
the  vicinity  of  the  edge.  The  Bessel  function  subroutine 
must  be  verified,  therefore,  for  a  range  of  argument 
0.001^  k^^  100  (Note:  Another  very  important  reason  for 
choosing  the  upper  bound  of  k^=  100  is  that  the  Bessel 
function  subroutine  requires  a  large  amount  of  computer 
memory  and  time  when  computing  large  orders  for  arguments 
greater  than  100  or  so.  Appendix  C  describes  the  basic 
theory  behind  the  Bessel  function  subroutine  used). 

liquations  (45)  and  (46)  for  the  wedge  involve  an 
infinite  summation.  Since  the  computer  can  only  evaluate 
finite  summations,  an  approximation  of  Hqs  (45)  and  (46) 
must  be  made  by  truncating  the  infinite  summation.  The 
truncated-off  terms  must  be  small  enough  so  as  to  be 
insignificant  if  all  of  them  were  added  to  the  finite 
Summation  approximation.  Recalling  that  the  Bessel 
function  of  concern  involves  an  order  having  an  integer 
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part,  N  ,  and  a  fractional  part.,  V  ,  it  is  necessary 

to  truncate  the  summation  at  some  maximum  integer  order, 

N  .  The  condition  chosen  on  N  is  such  that  for 

max  max 

any  integer  M  greater  than  N1  ,  the  absolute  value 

inn  x 

of  -Wkf)  must  be  less  than  10~10,  i.e. 

For  Jv*W(kfK 

where  0.001  100  ;  k^>  real 

0^  V  <  1  ;  V  real 

N  =  (0,1, 2, 3 . Nmax)  T75) 

then  for  any  integer  M  >  N 

'  °  max 

|'Wm  cke  >  |  < 10'10 

The  condition  for  |  f  )|  <  10"10  is  chosen 

tor  reasons  discussed  in  Appendix  B.  The  above  conditions, 

I 't  (75),  are  true  for  all  the  criteria  presented  in  this 
chapter. 

The  criteria  can  now  be  established  to  determine  the 
validity  of  the  Bessel  function  subroutine  for  the  range 
of  argument  and  order  given  in  Hq  (75). 


order  given  in  Eq  (75),  then  require 


XtnCMlsTP 

STD 

where  the  value  of  J  . ,  fk  p  )  calculated  by  the  sub- 

V+N  K  7 

routine  is  denoted  by  (k^>  anc*  t*ie  value 

of  Vn  (kf  )  given  by  accepted,  published  values  is 
denoted  by  Jv*M(k£)|STD  * 

The  value  of  10  relative  error  was  chosen  rather 
arbitrarily  since  it  is  difficult  to  predict  what  the 
computer  round-off  error  and  addition  of  absolute  error, 
due  to  the  summation  ir.  Eqs  (45)  and  (46),  will  come  out 
u>.  Round-off  error  is  defined  as  the  error  taking  into 
account  the  truncation  of  an  infinite  summation  and  the 
inherent  error  due  to  rounding  off  the  values  of  terms 
other  than  the  Bessel  function  in  Eqs  (45)  and  (46).  The 
round-off  error,  if  a  problem,  will  become  evident  in  the 
half-plane  check  in  the  next  section. 

As  a  check  to  see  how  the  absolute  error  of  the  Bessel 
function  subroutine  adds  up  in  a  summation,  use  is  made 
of  the  following  relations  (Ref  5:301,  361) 


<  10 


-5 


^  - 1 V 
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(76a) 


JF 


(76b) 


1«IM  +  2.Z  J,mM 


k  here 


£-U$*)  “nd  J  (Vf*1) 


are  commonly  known  as 


r he  Fresnel  cosine  and  sine  integrals,  respectively. 

•or  computational  reasons  discussed  earlier,  the  infinite 
sums  in  Eqs  (76a),  (76b),  and  (77)  must  be  truncated  to 
•\  ,  where  N  is  defined  as  in  Eq  (75).  Even  though 

I qs  (76a),  (76b),  and  (77)  involve  only  the  summation  of 
Bessel  functions  of  integer  order  or  of  integer  plus  one- 
ti.il  l~  order  while  Eqs  (45)  and  (46)  involve  the  summation 
of  Bessel  functions  of  any  fractional  order,  Eqs  (76a), 

1  ~(>b)  ,  and  (77)  provide  both  a  worst  case  and  best  case 
tor  x  ^.10  .  For  x  <10  ,  the  accuracy  will  depend 

on  computer  memory  requested  (see  Appendix  C) .  The  Bessel 
function  routine,  for  x  10  ,  provides  better  accuracy 

for  fractional  orders  approaching  one-half  than  for 
fractional  orders  approaching  zero  or  one  (see  Appendix  C) . 
As  a  result,  Eqs  (76a)  and  (76b)  provide  the  best  case  for 
checking  for  accuracy  and  Eq  (77)  provides  the  worst  case. 
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Criterion  (IV-2a).  For  the  range  of  argument  and 
order  given  in  Eq  (75),  then  require 


<t(^ke)  - 


<t  (vf«e) 


<  10 


-3 


where  £(7?V)  denotes  accepted,  published  values  of 

the  Fresnel  cosine  integral  and  the  summation  is  the 

truncated  version  of  Eq  (76a).  N  ^  in  this  case  is 

^  max 

equal  to  the  even  value  of  2N  where  the  summation  is 
t  runcated. 


Criterion  (IV-2b).  For  the  range  of  argument  and 
order  given  in  Eq  (75),  then  require 


N, 


ITmax 

t)  -  £  r 


sm 


M‘  f) 


<  Iff 


-3 


where  I  U/:E  kP 1  denotes  accepted,  published  values 
of  the  Fresnel  sine  integral  and  the  summation  is  the 
truncated  version  of  Eq  (76b).  N  ,  ,  in  this  case  is 
equal  to  the  odd  value  of  2N+1  where  the  summation  is 
t  runcated. 


Criterion  (IV-3).  For  the  ratine  of  argument  and 
order  given  in  Eq  (75),  then  require 


"maX 

-  [  X  M  +  2  I  JaN  (kf)l 

hi  *“  1 


<  to 
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where  the  expression  in  the  square  brackets  is  the 

truncated  version  of  Eq  (77).  N  in  this  case  is 

n  max 

equal  to  the  even  value  of  2N  where  the  summation  is 
r  runcated. 

As  in  Criterion  (IV-1),  the  relative  error  for  the 
most  recently  presented  criteria  was  chosen  rather 
arbitrarily  since  it  is  hard  to  predict  if  the  round-off 
error  will  be  a  problem.  The  next  section  (involving 
tiie  half-plane)  offers  a  better  indicator  on  how 
( in) significant  the  round-off  error  is. 


Results .  The  following  tables  arc  results  of  computer 
runs  to  test  the  criteria  presented  so  far  in  this  chapter, 
'nly  a  few  of  the  hundreds  of  values  tested  are  presented 
in  these  tables  (additional  tables  arc  in  Appendix  D) . 
li-.e  range  of  order  and  argument,  however,  is  limited  by 
the  non-availability  of  published  tables  to  compare 
against  (see  Appendix  D) . 

Table  (IV-1)  is  in  reference  to  Criterion  (IV-1) 

(two  additional  tables  are  in  Appendix  D) .  The  column 
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labeled  Relative  Error  gives  the  absolute  value  of  the 


error  given  in  Criterion  (IV-1).  As  evidenced  by 
1'able  (IV-1),  the  maximum  relative  error  between  the 
>ubroutine  value  and  the  accepted  value  is  less  than 
10  The  additional  tables  presented  in  Appendix  D 
also  have  maximum  relative  errors  less  than  10  All 

the  values  produced  by  the  subroutine  that  were  checked 
with  accepted  values  met  Criterion  (IV-1). 

Table  (IV-2a)  is  in  reference  to  Criterion  (IV-2a). 
lables  for  Criterion  (IV-2b)  and  Criterion  (IV-3)  are 
[■resented  in  Appendix  D.  The  column  labeled  Relative 
t  i ror  gives  the  absolute  value  of  the  error  given  in 
Criterion  (IV-2a).  As  evidenced  by  Table  (IV-2a),  the 
maximum  relative  error  between  the  subroutine  value  and 
the  accepted  value  is  less  than  10  ^ .  The  additional 
tables  presented  in  Appendix  D  also  have  maximum  relative 
errors  less  than  10  All  the  values  produced  by  the 
.-nbrout  ine  that  were  checked  with  accepted  values  met 
t  i  i ter i on  ( I V-  i )  . 


TABLE  (IV- la) 


Bessel  Function  Suhroi nine  Check 
Against  Published  Standards  for  Jy  +  (Sj(kf) 

Range:  0.01  ^  k  ^  100 

N  =  (0,1,2,...  ,100) 

V  -  1/2 


N 

Subroutine 

JN+i(ke  > 

Standard* 

Relative 

Error 

0 

.0797871199 

.0797871263 

8. 02(-8) 

1 

.  2659S8839S (-3) 

. 2659588606(-3) 

7.93(-8) 

2 

. 5319191988(-6) 

. 53191 924 1 1 ( -6) 

7.95(-8) 

12 

. 1009230707(-37) 

. 1009230776(-37) 

6. 84(-8) 

0 

.1783380682 

.1783580824 

7. 96(-8) 

1 

,2972796639(-2) 

.  2972 796874 (-2) 

7.91 (-8) 

2 

.2973009005(-4) 

.2973009240 (-4) 

7. 90 (-8) 

12 

. 55092972289(-29) 

.55092976611 (-29) 

6. 75 ( - 8) 

0 

.2518929203 

.2518929403 

7 . 94 ( - 8) 

1 

.8402033633(-2) 

.8402034300 (-2) 

3.96(-8) 

2 

. 1680887056(-3) 

. 1680887190(-3) 

7.97 (-8) 

12 

.3190882666(-25) 

.31 90882824 (-25) 

4 . 95 ( - 8) 

0 

.5409737469 

.5409737899 

7 . 95 ( - 8 ) 

1 

. 9170169233 ( - 1 ) 

.  91 701 099(H)  ( - 1 ) 

7. 93 (-8) 

2 

. 9236407085( -2) 

.9236407820(-2) 

7 . 96 ( - 8) 

12 

. 1734224713( -16) 

.  1 734224842 ( - 16) 

7.44(-8) 

0 

.6713966538 

.6713967071 

7. 94 (-8) 

1 

.2402978200 

.2402978391 

7 . 96 ( - 8) 

2 

. 4949680629 (-1) 

. 4949081 022( -  I) 

7. 94  ( -  8) 

12 

.9907033371 (-13) 

.99070341 82 (-13) 

8  19C-8) 

0 

-.1372637248 

-. 1372637358 

8.01(-8) 

10 

.1630073534 

.1630073664 

7. 98 (-8) 

0 

-.  2960583121 (-1) 

-. 29605831 89( - 1 ) 

2.30(-8) 

10 

-,8484972142(-1) 

- . 8484972092  f - 1 ) 

5. 89 (-9) 

100 

.  5749161075(-21) 

. 57491610 13 (-21) 

5. 57 (-9) 

0 

-.4040213234  (-1  ) 

-. 1010215272  1- 1 ) 

9.41 (-9) 

10 

-. 1561 124543(-2) 

-.1561 1 23855 (-2) 

4.41 (-7) 

100 

.  868 1 364822  C -  1 ) 

.  So815o4620( -  1 ) 

2. 33 (-8) 

*(IA;  Sf 

Not v :  Tlie  number  in  parenthesis  is  the  pov.vv  of  10  by  which  the 
preceding  number  must  be  multiplied. 
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TABLE  (IV -2a) 

Fresnel  Cosine  Integral  Check 
of  Subroutine  Versus  Published  Standards 


Range : 

0.0014137  ^  k 
N  =  (0,2,1, 

V  =  1/2 

^  ^  100.53096 

!) 

ke  -1  “2 

(^(U) -Subroutine  £(U) 

(see  Eq  (76a)  Standard* 

Relative 

Error 

.03 

.0014137 

.0300000 

.0300000 

0.0 

.04 

.0025133 

.0400000 

.0400000 

0.0 

.  l;6 

.0056549 

.0599997 

.0599998 

1 . 67  ( -6 ) 

.1)8 

.0100531 

.0799986 

.0799992 

7 . 50  (-6) 

0  .  1 

.015708 

.0999959 

.0999975 

1 . 6 2  (-  5  ) 

o . : 

.0628319 

.1998684 

.1999211 

2 . 64  ( -4 ) 

0.8 

.3926991 

.4923317 

.4923442 

2 . 54  (- 5  ) 

1.0 

1.5707963 

.7798923 

.7798934 

1 .41  (-6) 

1  .5 

3.5342917 

.4452353 

.4452612 

5 . 8 1  (  -  5  ) 

2.0 

6.2831853 

.4882271 

.4882534 

5 . 39 ( -  5 ) 

2.5 

9.817477 

.4574086 

.4574130 

5 . 62 (-6) 

3.0 

14.137167 

.60569] 1 

.  6  057208 

4 . 9  0  ( -  5  ) 

3  .  3 

19.242255 

.5325688 

.5325724 

6. 76 ( - 6 ) 

4  .0 

25.132741 

.498415 1 

.4984260 

2 . 1  3  ( -  5  ) 

4  .  5 

31.808626 

.5260180 

.5260259 

1 . 50  ( -  5 ) 

5 . 0 

39.269908 

.5636193 

.5636312 

2 . 1 0  (  -  5 ) 

6.0 

56.548668 

.4995256 

.4995315 

1  .  24 (- 5 ) 

o  .  5 

66.366145 

.4815931 

.4816035 

2 . 21  (-  5 ) 

7 . 0 

76.96902 

.5454600 

.5454671 

1 .35(-5) 

7.5 

88.357293 

.5160098 

.5160183 

1 . 7  0  ( -  5  ) 

8.0 

100.53096 

.4997858 

.4998022 

3 . 33 (- 5) 

*  (Ref 

5) 

Note  : 

The 

number  in  parenthesis  is 

the  power  of 

10  by  which 

the 

preceding  i 

number  must  be 

mu  1 1 i pi i ed . 

co 

J'J 


Conclusions .  Even  though  it  is  impossible  to  check 


for  every  possible  value  of  kp  ,  V  >  and  N  due  to 
the  infinite  number  of  values  kp  ,  V  ,  and  N  may  be 
(though  limited  in  range),  the  results  just  presented  for 
the  values  of  argument  and  order  chosen  all  satisfy  the 
criteria  of  this  section.  The  values  of  argument  and  order 
were  chosen  so  as  to  best  exhibit  the  behavior  (error-wise) 
of  the  subroutine  for  varied  ranges  of  argument  and  order 
(and  also  on  the  availability  of  published  standards).  The 
final  criterion  that  must  be  met  is  for  the  half-plane, 
presented  in  the  next  section. 

Criterion  and  Results  for  Verifying  the  Half-Plane 

The  purpose  of  this  section  is  to  verify  the  use  of  a 
truncated  version  of  Eqs  (52)  and  (52')  for  the  half-plane. 

By  the  similarity  of  Eqs  (52)  and  (52’)  to  Eqs  (45)  and 
(46)  for  the  wedge,  it  is  presumed  th3t  if  the  half-plane 
version  can  be  verified,  then  the  calculations  for  the  highly 
conducting  infinite  wedge,  presented  in  the  next  chapter, 
should  have  about  the  same  order  of  percent  error  as  the 
ha  1 f -plane . 

The  solutions  of  Born  and  Wolf  for  the  half-plane  (see 
Chapter  III)  are  used  as  the  standard  on  which  truncated 
versions  of  Eqs  (52)  and  (52’)  are  compared  for  relative 
error  determination. 


so 


Criterion  (IV-4) .  Use  is  made  of  Eqs  (55)  and  (55') 


from  Chapter  III  (Born  and  Wolf’s  classical  half -plane 
solution),  i.e. 


J  =  -H  a 

SBW  Z 

f 

where  H  =  P  _/« — ■*-  P  . ■ 

and  where  u  =  kf1  cos 

a  =  /JL  u 

*V  IT 

(55) 

\~v..  -  t  •  «... 

(55’) 

1 B IV 


ami  the  truncated  version  of  Eqs  (52)  and  (52') 
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where  the  range  of  argument  and  order  are  given  in  Eq  (75). 

N  in  this  case,  however,  is  equal  to  the  integer  portion 
of  m/2  where  the  summation  is  truncated.  It  is  then 
required  that 


be  less  than  0.01  for  excellent  results  and  less  than  0.05 
for  good  results.  Again,  as  before,  the  choices  for  0.01 
and  0.05  are  completely  arbitrary.  The  choice  normally 
depends  on  the  accuracy  needed  for  a  particular  application 
or  experiment. 

The  range  of  incident  radiation  angle,  (f  ,  is  chosen 
so  as  to  coincide  with  the  range  of  angles  of  incidence, 
measured  with  respect  to  the  surface  normal,  of  the  wedge 
presented  in  the  next  chapter. 


Results .  The  following  tables  are  results  of  computer 
runs  to  test  Criterion  (IV-4).  The  angle  of  incidence  is 
now  described  by  0;  where  0;  is  measured  with  respect  to 


the  surface  normal  (see  Figure  14).  Table  (IV-4a)  is  for 
0;  =  0°  (normal  incidence)  and  Table  (lV-4b)  is  for 
0j  =  45°  .  Additional  tables  for  0(  are  given  in 

Appendix  D.  Without  affecting  the  relative  error,  the 


quant i t  y 


is  normalized  to  equal  one. 
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TABLE  (IV- 4a) 


Half-Plane  Check  of 

Subroutine 

Versus  Bom  and  Wolf 

's  Solutions 

Angle  of  Incidence : 

e,  =  o' 

Distance  from  Edge : 

:  0. 001  100 

■** 

J 

s  s 

j.3* 

Relative 

k? 

Subroutine" 

Rom  and  Wolf* 

Error 

.001 

.3445668154(3) 

.3445668614(3) 

1 .63 (-7) 

.002 

.1780327007(3) 

.1780327274(3) 

1 . 50 (-7) 

.003 

.1217155868 (3) 

.1217156048(3) 

1.48 (-7) 

.004 

.9324531437(2) 

.9324532796(2) 

1.46(-7) 

.005 

.7600434579(2) 

.7600435669(2) 

1 . 44 ( -7 ) 

.006 

.6441669591 (2) 

.6441670484(2) 

1.39(-7) 

.007 

.5607918191(2) 

.5607918942(2) 

1.34 (-7) 

.008 

.4978420458(2) 

.4978421095(2) 

1.28(-7) 

.009 

.4485782759(2) 

.4485783299(2) 

1.20(-7) 

.01 

.4089399381(2) 

.4089399834(2) 

1.11  (-7) 

.05 

.1113554737(2) 

.1113547673(2) 

6.34 (-6) 

0.1 

.7004508694(1) 

.7004651448(1) 

6.24 (-5) 

0.5 

.3581315193(1) 

.3581192719(1) 

3.41 (-5) 

1  ,0 

.3425171544(1) 

.3125149724(1) 

6.37 (-6) 

5.0 

.4064911798(1) 

.4064741893(1) 

4.17 (-5) 

1  0 

.4032786222(1) 

.4032755312(1) 

7.66(-6) 

2  0 

.3988001132(1) 

.3987946918(1) 

1 . 36 ( - 5 ) 

50 

.4004309739(1) 

.4004269284(1) 

1.0K-5) 

4  0 

.3999589141(1) 

.3999588194(1) 

2 . 37 (-7 ) 

5  0 

.3998500232(1) 

.3998500575(1) 

8.57 (-8) 

60 

.4002128775(1) 

.4002128661(1) 

2.83  (-8) 

7  0 

.3998080504(1) 

.3998080486(1) 

4 . 54 (-9) 

8  0 

.4001249244(1) 

.4001249246(1) 

6.53(-10) 

90 

.3999562725(1) 

.3999562692(1) 

8 . 24 (-9) 

1  (>:■ 

.3999732579(1) 

.3999732559(1) 

5. 19(-9) 

Note:  Number  in  parenthesis  is  the  power  of  10  by  which  the  number 
must  be  multiplied. 


f  Normalized  for 
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TABLE  (IV -4b) 

Half-Plane  Check  of  Subroutine 
Versus  Born  and  Wolf's  Solutions 


Angle  of  Incidence:  8;  =  45* 

Distance  from  Edge:  0.001  ^  100 


1. 

j.y 

s  s 

Subrout ine1 

r 

Born  and  Wolff 

Relative 

Error 

.001 

.5672064751(3) 

.5672065603(3) 

1.50  (-7) 

.002 

.2886873766(3) 

.2886874211(3) 

1 . 54 ( - 7 ) 

.003 

.1951000809(3) 

.1951001107(3) 

1 . 53 (- 7) 

.004  • 

.1480164777(3) 

.1480165002(3) 

1 . 52 (-7) 

.005 

.1196182206(3) 

.1196182386(3) 

1.51(-7) 

.006 

.1005984697(3) 

.1005984848(3) 

1.50  (-7) 

.007 

.8695610556(2) 

.8695611858(2) 

1 . 50 (-7) 

.008 

.7668500398(2) 

.7668501540(2) 

1 . 49 (-7) 

.009 

.6866782750(2) 

.6866783768(2) 

1 . 48 (-7) 

.01 

.6223260214(2) 

.6223261132(2) 

1 . 48  (- 7) 

.05 

.1471591210(2) 

.1471591248(2) 

2 . 58 (-8) 

0.1 

.8339781086(1) 

.8339773058(1) 

9 . 62 (-7) 

0.5 

.2816598597(1) 

.2816001727(1) 

2 . 12 (-4) 

1.0 

.2068285593(1) 

.2068281692(1) 

2.37 (-6) 

5.0 

.1784461110(1) 

.1784323115(1) 

7. 73 (-5) 

10 

.2019055026(1) 

.2019047014(1) 

4 .46 (-b  J 

20 

.1995248743(1) 

.1995224135(1) 

1 . 23 ( - 5 ) 

30 

.1999774545(1) 

.1‘>99747669(1) 

1 .34  (-5) 

40 

.2002617145(1) 

.2002596108(1) 

1 . 05 ( - 5) 

50 

.1996340769(1) 

. 1996326280(1) 

7.25  (-6) 

60 

.2004134840(1) 

.2004125591(1) 

4.61 (-6) 

70 

.1995775856(1) 

.1995723727(1) 

2.61(-5) 

SO 

.2004112916(1) 

.2004083762(1) 

1.45 (-5) 

90 

.1996154167(1) 

.1996138054(1) 

8.07 (-6) 

100 

.2003493654(1) 

.2003429850(1) 

3. 18  (-5) 

t  Normalized  for  /  6a.  E  =1. 

Note:  The  number  in  parenthesis  is  the 

power  of  10  by  which  the  preceding 

number  must  be  multiplied. 
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Conclusions .  All  of  the  results  obtained  were  excellent 
(<0.01  relative  error  in  accordance  with  Criterion  (IV-4)). 

In  addition  to  the  tables  shown  in  this  chapter  and  Appendix  D, 
checks  were  made  for  varying  0;  iy  0.1  radians  (0^  Q-  ^  ^  ) 
with  the  same  results  of  "excellent.”  The  results  for  the 
infinite  half-volume,  given  in  Appendix  D,  also  came  out  as 
exce 1 lent . 

Since  Eqs  (45)  and  (46)  for  the  infinite  wedge  are  similar 
in  nature  to  Eqs  (52)  and  (52')  for  the  half-plane,  then  great 
promise  is  given  to  the  success  of  a  similar  truncation  working 
on  the  wedge  equations. 
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For  computational  reasons,  the  infinite  summations 
iven  in  Eqs  (46)  and  (49)  are  truncated  in  the  manner  of 
r.e  half-plane  case  verified  in  Chapter  17.  The  range  of 
ngles  of  the  incident  monochroinati c  plane  wave  is  chosen 
o4Q-,4  [id  (where  ©j  is  measured  a:;  show/,  in  Figure  15)  • 

:.is  range  is  chosen  because  it  provides  the  worst  case 
nighcr  power  density  loss)  for  the  top  wedge  surface  ( <f>~  ^  ) 
ear  the  vicinity  of  the  edge.  Thus,  the  range  of  incident 
..gles  for  the  incident  radiation  covers  normal  incidence 
9;  =0  )  to  grazing  incidence  (0^  =  90  ).  The  range  of  distance 
um  the  wedge  edge  times  the  propagation  constant  k  is 
-  Ml  ^  K  p  ^  100. 

:  nc-  following  graphs  yield  the  normalized  ( 

'.me -averaged  power  density  loss  (in  watts  per  square  meter) 

•  r  buth  TE  and  TM  polarization  of  the  incident  field.  I/.uiti- 
lica’icn  of  the  normalized  power  density  loss  by 
r  a  give*:.  O’  ,  £?  ,  and  E0  thus  yields  a  f:  rst  approximate 
t:  v.at*  per  square  meter  absorber  by  t:.e  wedge  r.e  a.  i:.v 
i  r. .  ty  wi  the  edge.  fr.e  cor.  :  i  t  i  >u  CT»U)€  is  assume 
■lisiied.  The  graphs  are  of  two  types:  (1)  normalizes  pa».er 
ensity  loss  versus  k^  for  various  values  of  0.  5  (2)  nor- 
alized  power  density  loss  versus  0f  for  various  values  of  \ 
itie  tables  give  the  power  density  loss  in  terms  of  an 
r  sorption  coefficient.  The  absorption  coefficient,  ^  ,  i_ 

:.e  ratio  cf'  absorbed  power  density  to  the  incident  power 
..sity  ( i rradiance ) j  i.e.  e =siy  1.  ,  where  ^  i s  given 


as  in  Eq  (18)  and  I0  is  the  irradiance  (in  watts  per  square 

«L 


meter)  and  equal  to  E0  , 
if  the  normalized  power  density  is  multiplied  by 


dince  the  graphs  yield  v  , 

60  Eq/o-^,. 


then 


n  (?u>.  (e.E  t/cr 

''  I.  "  VW  E 


TT 

-o 


<T S  V  JUt> 


(79) 


-nere  ^ p0rm  norrnalized  power  density  loss  (  as  in 

the  graphs).  The  results  are  for  both  T2  and  TM  polarization 
of  incident  field  wavelengths  of  3-8yUm  and  10.6yW.rn  and 
.C  n  e  conductivity  of  <T  =  5-8  x  10?  mhos/meter.  The  column 
’aoeled  0||max  gives  the  value  of  6{  such  that  (3  is  maxi¬ 
mum.  at  a  given  value  of  •  The  column  labeled  ^  (A) 
rives  the  distance  from  the  wedge  edge  in  angstroms.  This 
aiae  depends  on  the  incident  field  wavelength  (  k  =  2 TT /  X ) . 

The  graphs  and  tables  ar-e  now  presented.  Analysis  ui 
iesults  follows  the  graphs  and  tables. 
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v.n 


.1 

J 


Curves  from  top  to  coctx  ire  for: 


TABLE  (V-t) 


Absorption  Coefficient  for  Wedge 
Incident  Wavelength!  \=  3-8/im 
Skin  Depths  S  =  74.4  & 


ke 

*  /  a  \ 

TE 

TM 

0  IA; 

on 

mm 

& 

0.001 

6.0479 

0.0° 

0.01107 

45.0° 

0.94672 

0.005 

30.239 

0.0* 

0.01132 

45.0° 

0.32390 

0.01 

60.479 

0.0° 

0.01155 

45.0° 

0.20418 

0.05 

302.39 

0.0° 

0.01269 

36.9° 

0.07149 

0.1 

604.79 

o.o* 

0.01365 

25.4° 

0.04672 

0.5 

3023.9 

0.0° 

0.01696 

8.2® 

0.02283 

1.0 

6047 . 9 

0.0° 

0.01746 

0.0® 

0.02140 

5.0 

30239. 

16.8* 

O.O349I 

0.0® 

0.02498 

10.0 

60479. 

48.3° 

0.04583 

0.0° 

0.02479 

50.0 

302394. 

74.1° 

0.03539 

0.0° 

0.02460 

100.0 

604789. 

77-0° 

0.03388 

0.0® 

0.02460 

:-s»  G>  =  0.05110  at  0 

\  max  \ 

—  2^192 

.  A  (k^j  =  4.0)  and 

Oil  =22.5* 

ana* 

!,  =  0.94672  at  0^6. 0479  A  (kp  =  0.001)  and  OJ  =  45^ 
max  v  V  'nwx. 


For  ^  100  A  , 

(^max5  °*14523  at  100.0  A  (kp  :  O.ol6and  0<|^=  45. d 
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TABLE  ( V - 2 ) 


Absorption  Coefficient  for  Wedge 
Incident  Wavelength t  A6=10.6^>tm 
Skin  Depth «  £=124.3  A 


ke 

f(A) 

TE 

TM 

*L. 

Q  jrttax 

0.001 

16.870 

0.0° 

0.00663 

45. 0° 

0.56686 

0.005 

84.352 

0.0* 

O.OO678 

45. 0° 

0.19394 

0.01 

168.70 

O.C? 

0.00691 

45.0° 

0.12226 

0.05 

843.52 

O.rf 

O.OO76O 

36.9° 

0.04280 

0.1 

I687.O 

O.C f 

0.00818 

25.4° 

0.02798 

0.5 

8435.2 

0.0 

0.01016 

8.2° 

0.01367 

1.0 

16870. 

0.01046 

0.0° 

0.01281 

5.0 

84352. 

16.0 

0.02090 

0.0° 

0.01496 

10.0 

168704. 

48.3 

0.02744 

0.0° 

0.01 4b 5 

50.0 

843521. 

74.1° 

0.02119 

0.0* 

0.01473 

100.0 

1687042. 

77.0° 

0.02029 

0.0° 

0.01473 

i'lii  -  0.03060  at  D 

'  HIAM  V. 

£  67482.  A  (k p  =  4.0)  and 

Oil  =22.5° 

TMi  (T^  0.56686  at  ^  ft! 6 870  A”  ( k p  0.001)  and  0;| ^  =  45 . c| 
For  p  ^  100  A  , 

0*17273  at  p*100.  A  (kp  =  0.006)  and  G;|m=x  45. 


Several  observations  are  evident  from  the  graphs  just 
resented.  Figure  16  through  and  including  Figure  19  are 
are  for  the  incident  field  having  T K  polarization  (the  elec¬ 
tric  field  z-polarized).  Figure  20  through  and  including 
•igure  23  are  for  the  incident  field  having  TE  polarization 
i the  magnetic  field  z-polarized). 

Figure  16  plots  versus  k^>  for  the  range  of 

0.001  ^  kf^O.Ol.  The  five  curves  are  for  different  values 

of  6;  as  indicated  on  the  graph.  It  is  evident  that  the 
maximum  values  of  (V  for  each  9; =  constant  curve  occur 

No<*n 

closest  to  the  edge  at  kp  =  O.opi.  The  maximum  value  of 
curs  for  B-  =  45°  and  k^?  =  0.001.  It  is  also  evident  that 
the  minimum  values  of  (?Kiorrn occur  for  @,-=  0°  (normal  incidence) 

and  0;=  90°  (gracing  incidence).  Figure  18,  which  plots  (J' 

versus  Q;  for  different  values  of  k^)  (0.001^kp$  0.01), 

..hows  that  for  each  k^  =  constant  curve,  is  maximum 


fJorrr) 


Worn 

o 


t.  6;  =  4*5  and  minimum  at  @;  =  0  and  0,=  90  .  All  of  the 
curves  are  nonlinear.  This  is  to  be  expected  for  TM  polar- 
:  sation,  since  for  0;  =  45°  (wedge  angle  ~  2  o<  =  90°),  the 
surface  current  parallel  to  the  edge  will  increase  as  ^  J 
as  the  edge  is  approached  (Ref  3»  19). 

In  direct  contrast  to  the  TM  case  for  the  range  0.001 
^  kp^  0.01,  is  the  TE  case  for  the  same  range  of  k  ^  . 

Figure  23  plots  ^  ^versus  k^>  ( 0 . 001  ^  k p  $  0.01 ) .  For 

e. »  ""  ~ 


0°  and 


0;  =  22.5°»  C^efmincreases  linearly  with  k  . 
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For  0,=  6?.5°and  6;  =  90°  •  6^  decreases  linearly  with  kp>  . 
por  0,—  45°.  (5^  r  essentially  remains  constant  for  increas¬ 
ing  k  p  .  The  value  of  s  also  much  smaller  (by  around 

.'’O' to  80  times)  than  (^^for  the  TK  case,  as  in  Figure  16. 

This  is  to  be  expected,  however,  because  ^lor  ^in  the  TM  case 
is  proportional  to  (kp  )-^  while  (5j^  in  the  TE  case  is  not. 
Thus,  very  small  values  of  k  p  in  the  TK  case  (kp<<  1) 
dramatically  increases  (5^^.  As  evidenced  from  Figures 
20  and  23  for  the  range  0.001  ^kp^  0.01  and  O^0,-^9O  , 
ft  doesn't  change  as  dramatically  for  the  TE  case  as  for 
the  TK  case  in  Figures  16  and  18. 

Figure  1?  shows  that  6^ „,*,de creases  rapidly  from 
k  p  <  0.5  to  kp>0.5.  and  after  the  'overshoot'  present  around 
l<kp<5  damps  out,  remains  essentially  constant  for 

k p  >  10.  This  constant  value  for  (5^or^de creases  with  in¬ 
creasing  G;«  This  is  to  be  expected  since  the  reflectance, 
which  is  proportional  to  the  reflected  field  flux  over  the 

incident  field  flux„  increases  with  increasing  angle  of  incidence. 

Figures  21  and  22  plot  ^jor(n  versus  kp  (0.5$  k  p  $1  100) 
for  various  values  of  ©,■  for  TE  polarization.  In  direct 
contrast  to  the  TM  case  for  this  range  of  k  p  ,  (5^  does  not 
damp  out  rapidly  and  approach  a  constant.  Figure  21  shows 
that  for  0;=  0°  and  6-,=  45°,  (5^of  oscillates  like  a  sinusoidal 
function  times  a  decaying  exponential  function.  Figure  22 
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A  «  ,  6  O 

shows  similar  behavior  for  o;  =  22.5  .  67.5  ,  and  90  .  The 
maximum  values  of  (Y.  occur  for  &;  =  22.5°  and  k  p  ^5 .  The 

Norm  V 

r.  lnimum  values  of  (5  occur  for  9:  =  90*  .  Figures  21  and 

'•'Norm 

:  .for  the  TE  case  appear  to  oscillate  about  £j^orwpqual  to 
four  (except  for  0;  =  90°),  whereas  in  the  TM  case, 
decreases  with  ©;  for  k^>  1  .  This  is  to  be  expected  because 
the  reflectance  in  the  TM  case  changes  less  dramatically  for 
increasing  ©;  (until  ©,-  approaches  90°  )  than  for  the  TE 
ease  (Ref  9s  80-89). 

Figure  19  plots  ®  versus  &;  for  various  values  of 

Nerm 

^  (0.1$k<  100).  In  relation  to  Figure  18,  the  values  of 

0?.  is  greatly  reduced.  Also,  except  for  kasC.l,  the 

Ncfi'rt  l. 

carious  k  £>  -  constant  curves  tend  to  lose  the  bell  shape 
evident  in  the  curves  for  k^^O.Ol  in  Figure  18.  For  k£>  10, 
I'r.e  curves  are  nearly  the  same  (as  evidenced  by  the  k^  =10 
curve  being  almost  the  same  as  the  k ^  -  100  curve  in  Figure 
i.o).  This  Lrend  was  evident  from  the  data  obtained  for  any 
ilue  of  kp  in  this  range.  The  TE  case  is  again  much  different 
•  ir  ti.-i  same  range  of  k  ^  .  The  various  curves  for  0.1  ^ 

^^100  in  Figure  20  do  not  share  a  common  shape  as  the  curves 
in  Figure  19  for  the  TM  case  do.  The  highest  values  of  O^0 
appear  on  the  k  ^  =  10  curve. 

Figures  16  and  19  can  be  compared  to  similar  plots  for 
the  half-plane  (see  Figures  24  and  25  in  Appendix  D) .  From 
.  .ch  a  comparison,  it  is  evident  that  for  the  range  of  k p 


closest  to  the  edge  (0.001 ^  k  p  ^  0.01),  the  induced  current 
ir.  the  half-plane  is  considerably  higher  than  for  the  90° 
wedge.  This  is  expected  since  for  the  half-plane,  the  surface 
current  parallel  to  the  edge  m<  reasea  as  P  as  the  edge  >s 

approached.  For  the  wedge,  the  surface  current  increases  as 

-  V3 

P  (Ref  3s  19).  For  distances  further  from  the  edge,  however, 
^lorrn  for  the  wedge  are  on  the  order  o  the  half- 

nlane.  This  indicates  that  the  wedge  reduces  to  the  half- 
lane  problem  for  k^  far  enough  away  from  the  edge.  This  is 
.  •  ne  for  both  the  TE  and  TM  cases. 

Tables  (V-l)  and  (V-2)  specialize  the  results  for  mater- 
.  1  parameters  (T,  €  ,  and  JX  and  incident  wavelength,  . 

.no  skin  depth,  &  ,  is  determined  as  in  Eq  (12).  As  given 
in  Eq  (79),  the  tables  give  the  maximum  absorption  of  incident 
:  1  ux  ( watts/meter2)  for  a  given  distance  from  the  wedge  edge, 

P  ,  and  Q. !  .  Thus,  the  tables  are  in  answer  to  the  main 

V  ’  I  wax 

o:  jective  of  this  thesis  described  in  Chapter  I. 

..aoIc  (7-1)  is  for  the  incideiit  field  wavelength  equal 
3-  ytxm.  .,'ith  a  conductivity  CT  assumed  to  be  equal  to 
•  xl 0  mos/meter  ,  the  skin  depth  S  is  thus  equal  to 

A 

4  .•»  (see  Eq  (15)  )•  The  maximum  absorption  of  incident 
f:> Id  flux  is  for  a  maximum  value  of  ^  .  For  TE  incident 

c 

polarization,  (i  occurs  for  a  distance  D  —  24192.  A  from 

t:.e  e  ige  and  0,]  s  22. 5°  •  The  amount  of  time-averaged  power 
.  ..city  absorbed  is  thus  equal  to  approximately  five  percent 


of  the  incident  field  flux  (in  watts  per  square  meter).  Thus, 
for  an  incident  field  flux  (TE  polarized,  @;  =  22. 5°  #  ^0=  3»8yum) 
of  10°  watts/cm2,  the  copper  wedge  would  have  to  be  able  to 
absorb  5  x  10^  watts/cm2  ,  otherwise  it  might  be  damaged  a 
distance  2.4^rr.  away  from  the  edge.  For  an  incident  field 

of  wavelength  10.6^m  ,  Table  (V-2)  reveals  the  most  absorp¬ 
tion  takes  place  at  ^  6.7^/im  and  dt  =  22.5°  (  (3  —  0.03). 

These  percentages  are  a  percent  or  two  higher  than  for  copper 
of  any  geometrical  shape  for  the  3-8/<4m  to  10.6^<m  range 
of  incident  field  wavelength  (Ref  9*  88). 

The  TiY  case  is  much  different.  Tables  (V-l)  and  (V-2) 

-  how  Q>  as  high  as  0.94672.  This  means  that  almost  95  per¬ 
cent  of  the  incident  field  flux  is  absorbed  at  k  =  0.001 
oi.i  0;  =  0°  .  The  calculated  is  only  a  first  approxi¬ 
mation,  however.  It  is  based  on  the  assumption  that  the 
wedge  is  a  linear,  isotropic  medium.  For  k ^  =  0.001  ,  the 
.tuaee  from  the  edge  is  only  about  1 6 . 9  A  for  10.6  Uv„ 

a..l  6.0  A  for  =  3.8^<m.  Eincu  copper  has  a  radius  of 
a:. out  one  angstrom,  the  linear  assumption  surely  breaks  down 
at  k^)  -  0.001.  This  problem  is  not  evident  in  the  TE  case 

because  (i  occurs  for  a  distance  from  the  edge  thousands 
'  wax 

of  copper  atoms  away,  hence  accepting  linearity  of  the  medium 
a.;  a  valid  assumption.  It  is  now  assumed  that  the  linear 
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o  9 

assumption  be  valid  for  p  100  A  .  The  value  of  100  A  is 
arbitrary,  but  necessary  in  order  to  establish  a  valid  inter¬ 


pretation  of  the  results.  Thus,  for  p  -  100  A  ,  kp  must  be 
greater  than  0.006  for  \  -  10.6  ^ in  and  kp  must  be  greater 
than. 0165  for  X*  =  3*8  .  As  evidenced  from  Figures  16 

and  17,  the  maximum  ^  will  thus  occur  at  kp  =  0.006  ( p  - 
100  A*  )  at  10,6 yi/m  and  6;  =  45*  and  at  k ^  =  .0l65({)£r 
100  A  )  at  X0=  3«8  /tm  and  0,=  45®  For  X0=  3*8  ^m  , 


aoout  0.14^  which  is  much  more 


m  , 

than  (3  for  TE  polar- 


Aia  X 


i nation.  For  XD=  ,  however  ,  the  value  of 

implies  that  more  than  17  percent  of  the  incident  flux  is 
...  sorbed,  Fcr  very  large  values  of  incident  flux,  damage  to 


0 

the  wedge  surface  100  A  from  the  edge  would  be  more  likely 
for  the  TK  case  than  for  the  TE  case. 

The  analysis  given  above  for  the  graphs  and  tables  are 
ualy  a  first  approximation.  The  assumptions  made  included 
treating  the  wedge  medium  as  linear  and  isotropic,  using  the 
..  nrfuce  current  density  for  a  perfect  conductor  (C->°°),  and 
u.e  sr.in  depth  approximation.  It  is  also  assumed  that  all  of 
the  power  density  absorbed  is  convertible  to  heat  (in  watts 
per  square  meter) .  Another  assumption  made  is  that  the  in¬ 
finite  wedge  solution  can  be  used  to  approximate  a  finite 
wedge  near  the  vicinity  of  the  edge .  Of  note  here  is  that  for 
kp  large  (greater  that  10)  ,  the  effects  of  the  edge  on  the 


induced  surface  current  is  dramatically  reduced  for  either 
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or  TE  polarization. 


V I .  Conclusions  and  We  commend  at  ions 

Vhe  basic  theory  to  sol  ve  the  wave  equation  for  given 
■  !om>.  :.riec  has  been  investigated  in  this  report.  This  theory 

as  been  applied  to  solve  for  the  time -averaged  power  density 

0 

issioated  on  the  surface  of  t  highly  conductive  90  infinite 

edge  illuminated  by  a  monochromatic  plane  wave.  By  multi- 

iying  the  normalized  results  by  a  constant  depending  on  the 

ate  rial  parameters  of  the  wedge  (  0" ,  £  ,  jX  )  and  the  incident 

ive ’ s  frequency,  numerical  results  are  derived  for  a  conduct- 

n 

.  .rface  of  G  -  5-8  x  10  mhos/meter  and  incident 

ngths  of  \=  3 . 8yt/ rn  and  \  =  10.6^um  .  Though  the 

edge  used  in  the  calculation  is  assumed  infinite  in  extent, 

:ic  effect;.  of  the  edge  on  the  induced  surface  current  is  seen 
o  be  dramatically  reduced  just  a  few  wavelengths  av/ay  from 

i  •  0  t: !  lK  [Q  • 

. nis  :  tuuy  has  shown  that  for  the  electric  field  of  a 

.  .O'  . .  roma  :  i  c  plane  wave  polar  i  zed  parallel  to  the  sharp  t  .  •_ 

:  a  highly  conductive  90° wedge ,  a  high  percentage  (  17  per- 

O 

v.  nt  j  of  tiu;  incident  flux  can  be  absorbed  100  A  away  from 
ne  edge.  Though  many  assumptions  axe  made  to  arrive  at  this 
esule,  it  has  been  shown  that  it  is  a  good  first  approxima- 
ion.  Therefore,  this  study  provides  a  keener  insight  into 

:.e  amount  of  incident  flux  absorbed  by  a  highly  conductive 

0 

•  «  ..edge  near  the  vicinity  of  the  edge  for  both  TS  and  TM 
■  la:  izati'.n  and  various  angles,  of  incidence. 

H  '5 


iince  the  results  are  only  a  first  approximation,  it 
;  recommended  that  more  factors  be  taken  into  account  for 
>re  _  scura ie  results.  These  include  thermal  conductivity  of 
•  ,e.  wedge  o'  well  as  applying  Fourier  transform  techniques 
r  a  finite  wedge.  Fv'ore  study  into  finding  the  difference 
.  : . .itu:^  of  the  surface  currents  induced  on  the  surface 
1  -inite  conductive  wedge  as  opposed  to  one  of  infinite 
.  '  :  viy,  ,  is  also  recommended. 

...is  study  also  opens  the  door  to  the  need  for  better 
.  i  .  function  routines  to  calculate  high  fractional  order 
..  i  l.ish  a.,  ament  Bessel  functions  of  the  first  kind,  as 
:  e...  ed  by  the  shortcomings  of  the  Bessel  function  subroutine 
: sen.  red  in  Appendix  G. 

it  is  finally  recommended  that  more  study  be  given  to 
riving  for  the  eigenvalues  and  eigenfunctions  for  the  type 
■  /.  problem  encountered  with  the  cone  and  parabolic  cylinder. 

.■.is  i  important  because  of  the  increase  in  use  of  such  con.- 
..  ....  as  ..  .1X1  cij .  1  1  s.  .  ! .  jx  i  l  oj  j.  j  . 
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APPENDIX  A 


Definition  of  Certain  Functions 

The  purpose  of  this  Appendix  is  to  lend  support  to  various 
references  made  within  this  thesis  to  the  definitions  of  certain 
functions . 


Bessel  Functions  (Ref  1:460-464) 

Bessel's  equation  of  order  V  is 


x  4-  c*  -  v2)>’  =  0 

d*  dx 

Solutions  to  Eq  (A-l)  include 


Jv(x) 


y  (-1)  M _ 

h  ml  Cl)*"" 


f  C 

fe  m!(m-v)lUr“V 


(A-l) 


(A-2) 


(A- 3) 


where  Jr(x)  and  J  v(x)  are  Bessel  functions  of  the  first 
kind  of  order  v  and  -v  ,  respectively.  Bessel's  function 
of  the  second  kind  of  order  v  is  also  a  solution  to  Eq  (A-l) 
and  is  given  by 

N  (x)  -  Jy.  ^ -C..°.S  T •  v. ftL  (A 

v  Si«  (Vir) 
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Hankcl  functions  of  the  first  and  second  kind  of  order  v 


are  also  solutions 

to  Eq  (A- 1 ) 

and  defined  as 

H^Cx) 

u 

c* 

< 

X 

+ 

i  Ny(x) 

(A- 5) 

h'2’(x) 

=  JvCx)  - 

i  Ny(x) 

(A-6) 

respectively. 

Asymptotic  approximations  for  Eqs  (A-2)  and  (A-4)  for 
x— >0  are  given  by 


prov  ided. 
leading 


Re 


>  0 


For  large  arguments  (x— — >>«°) 


ti-rms  in  the  asymptotic  expansion  are  given  by 


CA-7) 

(A-8) 


,  the 


provided  J phase  (x)  |  <  fl*  .  Equations  (A-9)  and  (A-10)  into 
Eqs  (A-S)  and  (A-6)  thus  yield 
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(x)  - - ^ 

(Ti^ 

(A-ll) 

K-*oc 

(x)  - > 

ITT  •  v  -i  x 

J&- 1  e 

(A-12) 

X-*oo  . 

Associated  Legendre  Functions  (Ref  7:164,  193) 
The  associated  Legendre  equation  is 

(i  -  z2)^!L  -  2z -j1  +  [v(vn) 


-41u= 


0  (A-13) 


for  arbitrary  V  and  m  =  (0 ,1 ,2 ,3  , . . . )  .  Solutions  to 

Eq  i A-13)  include 


?»to  =  u2  -  i)"^  V 


(A- 14 ) 


where  P^(z)  is  called  the  associated  Legendre  function  of  the 
first  kind.  P (z)  is  a  solution  to  Eq  (A - 14)  for  m  =  0  and 
called  Legendre's  function  of  the  first  kind.  It  is  given 
as 


P  (z) 
v  v  * 


y  k)i  (l-  z\k 

vto  '  2  1 


sin  (vff) 

TT 


f  (k-i-v)Hk^!/i-z\k 

hi  wr  UJ 


(A- 15 ) 
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where  N  is  the  closest  integer  to  V  such  that  N  ^  V  . 


Hermite  Functions  (Ref  7:260,  284,  294) 


i’he  differential  equation 

4\  ~  2-z  T~  + 

J.r-  <*Z 


=  0 


(A 


has  Hermite  functions  for  solutions.  Hermite' s  function  of 
degree  is  given  by 


'!  (z) 

v 


£0*1  l  Li  A  •  z*) 

r(^)  * 12  *  J  ' 

fM2  §(¥*?}*) 


r(- 

where  ^(ot.’tfjz)  is  known  as  the  confluent  hypergeometric 
function.  It  is  given  by 

oo 


(A 


for  |z\  <  0O 


and  ?  (0,-1, -2, -3,...) 


(A 


For  large  values  of  argument  (z— ^  oo)  t  then 


h»(°I-5o°>  (-Zz^  (A'19) 

-i(v+l)iri  -v-1 

H-v-itu,TTr^  e  Uz)  tA-2°) 


For  V  equal  to  an  integer  n  ,  then  the  asymptotic  expressions 
are  approximated  by 

H  (iz )  - *  in(2z)n 

Z~»0O 

i 

H,n_1  (z)  - ►  (2z)  "n" 1 

Z— *  oo 

I 


(A~21) 

(A- 22) 
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APPENDIX  B 


Miscellaneous  Support  of  Thesis 

The  purpose  of  this  Appendix  is  to  lend  support  to 
various  references  made  within  this  thesis. 


Derivation  of  Source  Current  Density 


The  details  leading  to  the  derivation  of  Eq  (36)  is 
presented  in  this  section.  In  reference  to  Figure  7a  in 
Chapter  III,  the  current  filament  is  considered  to  be  an 
impulse  of  current  of  strength  I  (in  amps).  Figure  7b 
represents  a  vanishingly  small  area  surrounding  the  current 
filament;  small  enough  so  that  the  current  filament  can  be 
considered  as  a  current  sheet  (£=  p/  )  separating  the 
regions  for  and  p  <  ^ '  . 

Applying  Ampere’s  circuital  law  (Ref  1:33) 


(B-l) 


where  i*  is  the  total  current  on  the  surface  of  the  boundary 
enclosed  by  the  path  shown  in  Figure  l\x,  then  (Ref  1:34) 


a  .  r-*  -*■  n 

n  X  [H,  -  HnJ  -  J$ 


(B-2) 


is  obtained  for  the  surface  current  density.  It  is  noted 

i 

^  and  Hj  is  the  magnetic  field  on  the  Region  I 


^  A 
that  n  *  a 
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I 


! 

I 


i 

1 

t 

i 


side  of  the  surface  boundary  and  Hjj  is  the  magnetic  field 
on  the  Region  II  side  of  the  surface  boundary.  Both  Hy 
and  must  be  £>-polarized  because  Jg  is  z-polarized 

^  A 

and  n  =  .  Taking  the  Region  I  side  of  the  boundary  to 

be  located  at  p  =  p’-f  and  the  Region  II  side  of  the  boundary 
to  be  located  at  f  -  e-  ,  then  the  scalar  surface  boundary 
eouation  is  written  as 


Jz  =  H,Cf'+  )  -  H(j  (  p'-) 


(B-5) 


where  1 1 . (  ^40 
Region  T  and 
in  Region  II. 


represents  the  surface  magnetic  field  in 

represents  the  surface  magnetic  field 
Recalling  Eq  (2') 


>  — * 

Vx  E  =  - 1 oojuc  H  fB*4) 

and  since  H  at  the  surface  is  ©-polarized  and  E  is 
£ -polar i zed ,  Eq  (2')  is  then  written  in  scalar  form  as 


ill 

3f 

Substitution  of  Eq  (31)  for  E into  the  above  equation 
yields 


(B-5) 
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-  Yj  qv  (<  ^ 


coju 


>v 


-  —  £  av3^(ke'j  Hv2)(ke)sin[v^-^)]sin[^  *)]  j  f  >e' 


u>f 


(B-6) 


where  .1  (k#)  denotes 


r  C  2 ) 7 


3TV0*) 

-  ^  “  and  (k^)  denotes 

#  Applying  Eq  (B-3)  to  Eq  (B-6),  with 

op 

for  ^  and  for  ^  <  £>  and  then  setting 

^  ^  ,  yields  for  the  source 

j  *  2  '-°<)]s\nl'' (&-*)] 

z  lO/Xo  V 


•  [j^e'K2)V)  -  ut%e'Wk(1] 


(B-7) 


Using  the  relations  (Ref  1:402) 


hJ2)(x)  =  Jv(x)  -  i  Ny(x) 
Hv2)/(x)  =  Jv(x)  "  1  Nv(x) 


(B-8) 


where  the  primes  denote  the  first  derivative  with  respect 
to  the  variable  x  and  N  (x)  is  the  Bessel  function  of  the 
second  kind  of  order  v  ,  and  then  by  applying  the  h’ronskian 
of  Bessel's  equation  given  by  (Ref  1:463) 
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tvoon;m  &  ( 

implies  that 

X  (*(')£%)- nfMx'M 

-  ivOcf'jivVi4- 1  V) 
=  -  i  Drv  (kf  1  K0<t ')  -  Nv  (kf 0  X  (kf ')] 


Substituting  this  result  into  the  expression  for  Jz  yields 

3 i=  “wr  (B  U) 

Lcontov.ch  Boundary  Condition  (Ref  3:18) 

As  long  as  the  radius  of  curvature  of  a  good  conductor 
is  much  larger  than  the  skin  depth,  the  conductor  boundary 
may  be  described  by  an  impedance  boundary  condition,  where 
f  of  hq  (17)  is  the  surface  impedance.  Since  the  wedge 

os 

has  a  flat  top  surface  and  the  point  of  observation  isn't 
on  the  tdge  singularity,  the  above  stated  l.eontovich  Boundary 
Condition  is  applicable. 
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Parabolic  Cylinder 


The  purpose  of  this  section  is  to  derive  Eqs  (62), 
( o 3 ) ,  and  (64).  Recalling  Eqs  (60)  and  (61), 


!(u**v*) 


ft  .  J*? 


m 

(3u 


2  + 


ft 

hr1 


f  kV=0 


jf'(ujVjz')  =  U(«)  VWZ(z) 


(60) 


(61) 


then  the  method  of  Separation  of  Variables  can  be  used. 
Substitution  of  Eq  (61)  into  Eq  (60),  dividing  through  by 
'f' t  0  ,  and  rearranging  yields 

i  fU  .  1  fV 


7  n!z 

V  avV  k  “~2  “s? 


(B-12) 


The  right  side  of  Eq  (B-12)  is  a  function  in  z  only,  while 

the  left  side  is  a  function  in  u  and  v  .  Therefore,  in 

order  for  the  above  equation  to  be  true  for  any  u,  v  , 

and  z  ,  both  sides  must  be  equal  to  a  constant.  Denoting 

2 

this  constant  by  kz  implies 


i  £Z  _  i 

"  Z  '  *z 


(B- 13) 


Substitution  of  Eq  (B-13)  into  the  previous  equation,  multiplying 
2  2  2 

through  by  c  (u  +  v  )  ,  and  rearranging  yields 
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The  right  side  of  Eq  (B-14)  is  a  function  in  v  only,  and 
the  left  side  is  a  function  in  u  only.  Therefore,  both 
sides  must  be  equal  to  a  constant  in  order  for  the  above 
equation  to  be  true  for  any  u  and  v  .  Denoting  this 
constant  by  1  ,  the  left  side  of  the  above  equation 

becomes 


(B-1S) 


and  the  right  side  becomes 

+  CB-w 

Since  Eqs  (B-13),  ( B - 1 5 )  ,  and  (R-]6)  are  functions  in 
only  one  variable,  the  partial  derivatives  are  written  as 
orJinar\  derivatives.  Also,  by  multiplying  Eqs  (B-13), 

(B-15),  and  (B-16)  through  by  Z  ,  U  ,  and  V  ,  respectively, 
and  rearranging  yields 


-o 

£V  +  [(c  K)V-(2/t*lj]Vi 

Cl  ^ 

**  /-s 


(B-17) 


(B-18) 


ipf.  f[(ckp)V+  (Zyu+iflv--  0 


2 

where  kn  is  defined  by 


k2  “  kD  » 


(B-19) 


A  substitution  of  variables  is  now  performed  so  as  to 
put  Eqs  (B-19)  and  (B-20)  into  the  form 


cUz 


[i  -  4 ‘]f(t) 


(B-20) 


which  has  parabolic  cylinder  functions  as  solutions  for 
and  where  the  variable  £  and  parameter  V  are,  in  general, 
real  or  complex  (Ref  7:284). 

Examining  Eq  (B-17)  first,  the  following  substitution 
of  variables  are  made: 


Then 


Q>  =  -Jl  ik»c'  “ 

(i2=  2  ikpc  U2 


Using  the  chain  rule,  then  7 — 

clu 


J.U. 


(B-21) 


ike 


MaKing  the  above  substitutions 


£  i  kDc 


into  Eq  (B-17)  yields 

[&p*-(wiu*o 


Dividing  the  above  equation  by  2ik^c  (^  0)  and  defining 


(B-22) 


results  in 


Jlfi 


(B-23) 


In  a  similar  fashion,  Eq  (B-18)  becomes 


where  the  substitution  of  variable 


(B-24) 


Y-  -^2  'K  c  V 

was  made  along  with  Eq  (B-22). 
Using  the  substitutions 


(B- 25 ) 


(3  -  -/i’u 

)T  =  V2’  * 
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(B- 26 ) 


(63) 


r 


into  Eqs  (B-23)  and  (B-24)  thus  results  in 

dlM  +  (zvf  i  -u2;u  =  o 
+  [-(z^l)-v*]V*0 

Criteria  Choice  o£  N 
_ _ max 

The  purpose  of  this  section  is  to  justify  the  choice  made 
for  N  in  Eq  (75). 

By  examining  Bessel  function  tables,  it  is  seen  that  the 
leading  terms  of  ar.  array  of  Jy+n(x)  (0  ^  1  ; 

_  ? 

0,1,2,...)  have  an  absolute  value  greater  than  10 
It  is  also  known  that  Jv+n(x)  is  a  decreasing  function 
for  V  +  f\  > x  (Ref  5:365,  368).  Thus,  by  choosing  to 
truncate  Eq  (45)  where  the  larger  ordered  Bessel  functions 
have  an  absolute  value  less  than  10  results  in  the 
truncated-cff  terms  (summed  together)  being  less  than  one 
percent  of  the  leading  terms. 
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APPENDIX  C 


Theory  of  Bessel  Function  Subroutine 


The  purpose  of  this  Appendix  is  to  briefly  describe 
tie  theory  behind  the  Bessel  function  subroutine  used  in 
this  thesis. 

The  subroutine  to  calculate  J  (x)  where  1 

y+n 

and  n  =  (0,1,2,...)  (y,  x  real)  is  a  UNIVAC  routine 

(Ref  10).  It  is  divided  into  three  regions  of  the  argument; 
(1 )  x  <  10  ,  (2)  10  £  x£50  ,  and  (3)  x  >  50  . 

(1)  x  <  10  (Ref  11:102-108) .  For  a  function  F  that 

obeys 


*  F„,<«)  *  ¥  F-"l 

then  let 

F.  (*)  =  *  +(3IW*) 

vf  n  ' 


for  M»  N  .  Defining 

F/+MH  ^  ~  ® 

F*hW  =  « 


(a  is  any  real  constant)  (C-2) 


it  can  be  shown  that  for 
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*  f;fh  60  -  °<  ^n(K) 


(c- 3) 


r 


n  >> 


X 

I 


and  using  the  addition  theorem 


L/Z 

2  1  IUf»M  =  « 

m=o 


i.t  re 


arra^j 

(vf  2rn)(v  +  M~i) 
•vi  (vf  am-2) 


(C-4) 


i’l.st  i  t  ut  ion  of  Eq  (C-2)  into  Eq  (C-l)  yields  an  array  of  F 

.ilaes  irom  F  (x)  to  F  tJj,(x)  .  It  is  desired  to  have 
v  V  ♦w'*  A 

M  >>  -g  so  ,:.q  (C-l)  can  he  used  to  recur  backward  without 
>>s  oi  accuracy.  Then  Eq  (C-4)  is  used  to  find  o(  >  and 
C*  in  l  q  (C.-3)  yields  ^y  +  n(x)  •  The  accuracy  will  depend 
vi..  mow  1  irge  M  is  chosen. 


(21  10^  x  ^  50  (Ref  12:18-26).  .J^(x)  can  be  described 


..■re  t  -  1/x  and 


104 


M*)  _ 

t 

t  \jk,  *  (  vf  7)  f\]  =  1  +  ^  +  ^  *■  *■  " 


tC-6) 


f  r  he  routine  truncates  at  t1^  ).  For  \vj  £  1  ,  Eqs  (C-5) 

and  (C-o)  yield  eight  figure  accuracy  in  Jv(x)  for  x  >  8 

Equation  (C-6)  is  also  most  rapidly  convergent  for  V  close 

:.c  one-half  (i.e.,  most  accurate  for  y-»  \  and  least 

£ 

accurate  for  y-*  0  and  y  1  ).  Eq  ( C-3 )  is  then  used 

to  recur  down,  thus  producing  an  array  of  Fv+n(x) 

Equation  (C-3)  is  then  used  to  find  o(  from  J^(x) 

vibst itution  of  0^  into  the  F^+n(x)  array  then  produces  a 

Jy+n(x)  array  (n  =  0,1,2 . N)  . 

(3;  x  >  50.  The  same  method  (phase-amplitude)  is  used 
for  x  >  SO  as  in  10 <  x  s£  50  .  However,  since  computer 

time  and  memory  are  limited  (especially  for  students),  it  is 
somewhat  more  difficult  to  get  M»  x/2  necessary  for 
backward  recursion.  Hence,  Jv(x)  and  J^  +  j(x)  calculated 
i.sing  Eqs  (C-5)  and  (C-6)  are  used  in  Eq  (C - 1 )  to  recur 
forward.  This  method  gives  poor  accuracy  for  large  orders 
because  for  ^(»x/2  ,  JyU)  —  (x/2  )^  . 


Modifications .  The  range  x  >  50  proved  to  be  unaccept¬ 
able  for  producing  accurate  values  of  Jy+n(x)  for  large 
n  O.  >  x)  .  Therefore,  it  was  decided  to  recur  backward 
for  the  extended  range  10<x  <100  .  This  proved  to  be 

very  successful.  Another  modification  attempted  involved 
calculating  a^  and  in  Eq  (C-6).  This  had  very  little 

effect  on  improving  the  routine. 

This  Bessel  function  subroutine  is  available  by  contacting 
Major  Glenn  R.  Doughty 
Air  Force  Institute  of  Technology 
Department  of  Physics 
Wright-Patterson  AFB  OH  45433 
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APPENDIX  D 

Additional  Numerical  Results 

Thi.~  appendix  provides  space  for  additional  figures  not 
placed  in  the  main  text.  The  non-availability  of  Bessel 
function  tables  for  large  order  and  argument  limited  the 
ranges  of  argument  and  order  that  could  be  verified.  For 
o ample,  Reference  13  only  provided  tables  for  fractional 
orders  limited  to  £  1/2  ,  £  1/3  ,  1  2/3  ,  and  +  3/4 
ii ! id  onl\  for  N  *  0  .  Reference  14,  on  the  other  hand, 
vw.s  limited  to  N  f  135  and  argument  <  100  ,  with 


i.  actional  order  equal  to  zero. 


TABLE  (IV -lb) 


Bessel  Function  Subroutine  Check 


Against  Published  Standards  for  J 

Range:  0.001  £  100 

N  -  (0,1,2,. . . ,100) 

V  =  0. 0 


v+ 


NCkf) 


k  c 
< 

N 

jN(kf) 

W 

Relative 

Subroutine 

Standard* 

Error 

.001 

0 

.9999997500 

.9999997500 

0.0 

.001 

1 

.0004999999 

.0004999999 

0.0 

n 

L 

.012499990(-5) 

.0124 99990  (-5) 

0.0 

3 

.0000020833 (-5) 

.0000020833  (-5) 

0.0 

0 

.9999750002 

.9999750002 

0.0 

.01 

l 

.0049999375 

.0049999375 

0.0 

2 

.0000124999 

.0000124999 

0.0 

3 

.0020833203(-5) 

,0020833203(-5) 

0.0 

0 

.7651976866 

.7651976866 

0.0 

1.0 

5 

.2497577302 (-3) 

.2497577302  (-3) 

0.0 

10 

.2630615124 (-9) 

.2630615124 (-9) 

0.0 

20 

.3873503009(-24) 

. 3873503009 (-24) 

0.0 

0 

-.1775967713 

-.1775967713 

0.0 

c  n 

5 

.2611405461 

.2611405461 

0.0 

o  ♦ 

10 

.1467802647 (-2) 

.1467802647 (-2) 

0.0 

20 

,2770330052(-10) 

.2770330052 (-10) 

0.0 

0 

-.2459357645 

-.2459357645 

0.0 

10 

10 

.2074861066 

.2074861066 

0.0 

20 

.1151336925 (-4) 

.11 5133692  5 ( -  4 ) 

0.0 

30 

.1551096078 (-11) 

.1551096078 (-11) 

0.0 

0 

.5581232812 (-1) 

.5581232767  (-1) 

8 . 06 (-9) 

50 

20 

-.1167043523 

-.1167043528 

4.28(-9) 

50 

.1214090214 

.1214090219 

4 . 12  (-9) 

100 

.1115927365(-20) 

.1115927368 (-20) 

2.69C-9) 

0 

,1998585061(-1) 

.1 998585030  (-1) 

1 . 55  (-8) 

100 

20 

.6221745833(-1) 

.6221745850(-1) 

2.73(-9) 

50 

3869834004 (-1) 

- .3869833973(-l) 

8.01(-9) 

- snrrm 

100 

.9636667398 (-1) 

n1  a1 /w  nsvnxv 

.9636667330  (-1) 

7.06(-9) 

- 

k p  =  b.OOl  and  0.01  (Ref  14);  For  ail  other  kf>  (Ref  5) 


Note:  The  nunber  in  parenthesis  is  the  power  of  10  by  which  the 
preceding  number  must  be  multiplied. 
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TABLE 

(IV-lc) 

Bessel  Function 

Subroutine  Check 

Against  Published  Standards  for  d^+N(kp) 

Range:  0.001  ^kp  ^100 

N  = 

0 

V  = 

1/4 

CL. 

-A 

J1/4W 

Subroutine 

J1/4(kp) 

Standard* 

Relative 

Error 

.001 

.1649761767 

.1649762131 

2.21C-7) 

.005 

.2466957406 

.2466957951 

2.21  (-7) 

0.01 

.2933679292 

.2933679941 

2.21C-7) 

0.05 

.4384768320 

.4384769287 

2.21(-7) 

0.1 

.5206577608 

.5206578756 

2.21(-7) 

0.5 

.7416564065 

.7416565702 

2.2l(-7) 

1.0 

.7522311674 

.7522313333 

2.21(-7) 

5.0 

-.2809720038 

-.2809720658 

2.21(-7) 

10 

-.2063937413 

-.2063937869 

2.21C-7) 

20 

.1782983389 

.1782983385 

2.24(-9) 

30 

-.1246044299 

-.1246044300 

8.03(-10) 

40 

. 5491175188 (-1) 

.5491175240C-1) 

9.47  (-9) 

50 

. 1410606328  (-1) 

. 141 0606269 (-1) 

4.18(-8) 

tU) 

-. 6642673498  (-1) 

- .6642673439  ( -1) 

8.88(-9) 

70 

.9125954031 (-1) 

. 9125953997 (-1) 

3.73(-9) 

80 

-. 8570889767 (-1) 

- . 8  57  0889765 ( - 1) 

2.33(-10) 

90 

.5511011408(-1) 

.5511 01 1439(-1) 

4. 72(-9) 

100 

-.1107092717  (-1) 

- ,1107092760(-1) 

3.88(-8) 

*  (Ref  13) 


Note.  The  number  in  parenthesis  is  the  power  of  10  by  which  the 
preceding  number  must  be  multiplied. 


TABLE  (IV- 2b) 

Fresnel  Sine  Integral  Check 
of  Subroutine  Versus  Published  Standards 


Range:  0.0014137  ^  k^>  ^  100.53096 

N  =  (1,3,S,. . . ,Nmax) 

V  =  1/2 


u 

kf  =  £u2 

(U) -Subroutine 
(see  Eq(76b)) 

CU) 

Standard* 

Relative 

Error 

.03 

.0014137 

.0000141 

.0000141 

0.0 

.04 

.0025133 

.0000335 

.0000335 

0.0 

.06 

.0056549 

.0001131 

.0001131 

0.0 

.09 

.0100531 

.0002681 

.0002681 

0.0 

0.1 

.015709 

.0005236 

.0005236 

0.0 

0.2 

.0628319 

.0041871 

.0041876 

1 . 19  ( -4) 

O.S 

.3926991 

.0647320 

.0647324 

6 . 18  (-6) 

1.0 

1.5707963 

.4382590 

.4382591 

2 . 28  (-8) 

1.5 

3.5342917 

,6975009 

.6975050 

5. 88 (-6) 

2.0 

6.2831853 

.3434101 

.3434157 

1 . 62  ( - 5 ) 

2  .  5 

9.817477 

.6191806 

.6191818 

1 . 62  ( - 6) 

3.0 

14.137167 

.4963041 

.4963130 

1 . 79  (-5) 

3 . 5 

19.242255 

.4152469 

.4152480 

2 . 65  ( -6) 

4 . 0 

25.132741 

.4205120 

.4205158 

9.04  (-6) 

4  .  S 

31.808626 

.4342700 

.4342730 

6 . 91  (-6 ) 

5.0 

39.269908 

.4991865 

.4991914 

5. 82  (-6) 

6 . 0 

56.548668 

.4469581 

.4469608 

6 . 04  ( -  6 ) 

6.5 

66.366145 

.5453716 

.5453764 

8.80  (- 6 ) 

7.0 

76.96902 

.4997014 

.4997048 

6 . 80  ( -6) 

7.5 

88.357293 

.4606970 

.4607012 

9.12  (-6) 

8.0 

100.53096 

.4602057 

.4602142 

1 . 91  (-5) 

*  (Ref  5) 

Note:  The  number  in  parenthesis  is  the  power  of  10  by  which 

the  preceding  number  must  be  multiplied. 
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TABLE  (IV- 3) 

Bessel  Summation  Check  of  Routine  for  Eq  (77) 


Range:  0.001  k^^lOO 

N  =  (0.1,2,. ...N  ) 

V  =  0 


kf 

Bessel  Sum 

A/in<u 

iOiJo(kf).2  S  J2N(k  ) 

Relative 

Error 

.001 

.99999975 

2 . 50  ( - 10) 

.002 

1.00000000 

0.0 

.003 

1.00000000 

o 

o 

.004 

1.00000000 

0.0 

.  005 

1.00000000 

o 

o 

.006 

1.00000000 

0.0 

.007 

1.00000000 

o 

o 

.  008 

1.00000000 

o 

o 

.009 

1.00000000 

o 

o 

.01 

1 . 00000000 

0.0 

.05 

.99999997 

3 . 00  (- 10) 

0.1 

1.00000000 

0.0 

0.5 

1  .00000000 

0.0 

1.0 

1 . oooooooo 

0.0 

5.0 

1.00000000 

0.0 

10 

1 . oooooooo 

0.0 

20 

1.00000000 

0.0 

30 

1.00000001 

1.00(-10) 

40 

1.00000000 

0.0 

50 

1.00000000 

o 

o 

60 

1.00000000 

0.  0 

70 

1.00000000 

0.0 

80 

1 . oooooooo 

0.0 

90 

1.00000000 

0.0 

100 

l . oooooooo 

0.0 

Note:  The  number  in  parenthesis  is  the  power  of  10  by  which  the  preceding 
number  must  be  multiplied. 


TABLE  (IV -4c) 

Half-Plane  Check  of  Subroutine 
Versus  Born  and  Wolf's  Solutions 


£ 

Angle  of  Incidence:  0;=  22.5 

Distance  from  Edge:  O.OOi  <  £  100 


ke 

J  *J* 
s  s 

Subroutine 

**  -**  + 

J  *J 
s  s 

Rom  and  Wolf 

Relative 

Error 

.001 

.4684029211(3) 

.4684029950(3) 

1.57 (-7) 

.002 

.2403209447(3) 

.2403209813(3) 

1.52(-7) 

.003 

.1634162387(3) 

.1634162633(3) 

1.50 (-7) 

.004 

.1246239813(3) 

.1246239998(3) 

1.49  (-7) 

.005 

.1011752171(3) 

.1011752320(3) 

1.47  (-7) 

.006 

.8544027553(2) 

.8544028804(2) 

1.46(-7) 

.007 

.7413464357(2) 

.7413465434(2) 

1.45  (-7) 

.008 

.6560952897 (2) 

.6560953843(2) 

1 . 44  (-7 ) 

.  009 

.5894561983(2) 

.5894562825(2) 

1.43(-7) 

.01 

.5358948389(2) 

.5358949148(2) 

1.42  f-  7) 

.05 

.1366385829(2) 

.1366384225(2) 

1.17  (-6) 

0.1 

.8189107747(1) 

.8189002466(1) 

1 . 29 (-5) 

0.5 

.3467604171(1) 

.3467593408(1) 

3.10(-6) 

1.0 

.2963554436(1) 

.2963273169(1) 

9.49  (-5) 

5.0 

.3469644246(1) 

.3469617950(1) 

7 . 58 (-6) 

10 

.3388221165(1) 

.3388136336(1) 

2. 50 (-5) 

20 

.3404785885(1) 

.3404684726(1) 

2 . 94 ( -5 ) 

30 

.3409823380(1) 

.3409763835(1) 

1.75(-5) 

40 

.3412019465(1) 

.3411990666(1) 

8.44 (-6) 

50 

.3413174576(1) 

.3413084812(1) 

2.63(-5) 

60 

.3413853061 (1) 

.3413817358(1) 

1 . 05 ( - 5) 

70 

.3414279933(1) 

.3414280489(1) 

1.63 (-7) 

80 

.3414559925(1) 

.3414560236(1) 

9.11 (-8) 

90 

.3413747647(1) 

.3414747812(1) 

4.87 (-8) 

100 

.3414873892(1) 

.3414873964(1) 

2.13 (-8) 

*  Normalized  for  _/i2?  =  1. 

^ o  O 

Note:  The  number  in  parenthesis  is  the  power  of  10  by  which  the 
preceding  number  must  be  multiplied. 
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TABLE  (IV-4d) 

Half-Plane  Check  of  Subroutine 
Versus  Born  and  Wolf's  Solutions 


Q 

Angle  of  Incidence:  6;  =  67.5 

Distance  from  Edge:  0.0f)l  <  100 


ke 

js-j* 

Subroutine 

t-J* 

s  s 

Born  and  Wolf 

Relative 

Error 

.  001 

.6259299676(3) 

.6259300619(3) 

1.51  (-7) 

.002 

.3158128663(3) 

.3158129157(3) 

1 . 56  (-7) 

.003 

.2120104401(3) 

.2120104731(3) 

1.56(-7) 

.004 

.1599423065(3) 

.159942331(3) 

1 . 55  (-7) 

.cos 

.1286161098(3) 

.1286161297(3) 

1.55(-7) 

.006 

.1076814883(3) 

.1076815049(3) 

1.54(-7) 

.007 

.9269540169(2) 

.9269541590(2) 

1.53C-7) 

.008 

.8143313033  (7) 

.8143314274(2) 

1.52C-7) 

.009 

.7265710850(2) 

.7265710602(2) 

1.5l(-7) 

.01 

.6562384447(2) 

.6562385430(2) 

1.50(-7) 

.05 

.1429570778(2) 

.1429570995(2) 

1.51 (-6) 

0.1 

.7620296496(1) 

.7620297497(1) 

1.31(-7) 

0.5 

.1995471355(1) 

.1995466379(1) 

2.49(-6) 

1.0 

.1219722187(1) 

.1219690885(1) 

2.57  (-5) 

5.0 

.5581411230 

.5581361479 

8.91  (-6) 

10 

.5005723246 

.5005719124 

1.02 (-6) 

20 

.5259130666 

.5258607850 

9.94 (-5) 

30 

.5671748589 

.5671596962 

2 . 67 (-5) 

40 

.5943320126 

.5942380397 

6.68(-6) 

50 

.6039451152 

. 6039073894 

6.25(-5) 

60 

.6004248186 

.6004162077 

1.43(-5) 

70 

.5909307243 

.5909287139 

3.40(-6) 

80 

.5820611890 

.5820490309 

2.09(-5) 

90 

.5777673975 

5.777645033 

5.01(-6) 

100 

.5787210245 

.5787067379 

2.47  (-5) 

"^Normalized  for  E_  =*  1 . 

W  jUe  O 

Note:  The  number  in  parenthesis  is  the  power  of  10  by  which  the 
preceding  number  must  be  multiplied. 
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TABLE  (IV -4e) 

Half-Plane  Check  o f  Su brout ine 
Versus  Born  and  Wolf’s  Solutions 


Angle  of  Incidence:  =  90 

Distance  from  Edge:  0.001  k  100 


k( 

J  -J* 
s  s 

Subroutine 

J  -J* 
s  s 

Born  and  Wolf 

Relative 

Error 

.001 

.6366196711(3) 

.6366197723(3) 

1 ,59(-7) 

.002 

.3183098356(3) 

.3183098861(3) 

l.S9(-7) 

.003 

.2122065570(3) 

.2122065908(3) 

1.59(-7) 

.004 

.1591549178 (3) 

.1591549431(3) 

1.59(-7) 

.005 

.1273239342(3) 

.1273239545(3) 

1.59(-7) 

.006 

.1061032786(3) 

.1061032954(3) 

1.59(-7) 

.  007 

.9094566739(2) 

.9094568175(2) 

1.59(-7) 

.  008 

.7957745888(2) 

.7957747153(2) 

1.59(-7) 

.009 

.7073551900(2) 

.7073553025(2) 

1.59(-7) 

.01 

.6366196710(2) 

.6366197722(2) 

1.59(-7) 

.05 

.1273239342(2) 

.1273239544(2) 

1 . 59  (-7) 

0.1 

.6366196711  (1) 

.6366197721(1) 

1 . 62 (-7 ) 

0.5 

.1273239338  (1) 

.1273239544(1) 

1 . 57 ( -7 ) 

1.0 

.6366195720 

.6366197718 

1.61(-7) 

5.0 

.1273239336 

.1273239542 

1.57 (-7) 

10 

.6366196710 (-1) 

.6366197709 (-1) 

1.69(-8) 

20 

.3183098905 (-1) 

.3183098851 (-1) 

5.47 (-9) 

30 

.2122065911 (-1) 

.2122065899 (-1) 

1 . 57 (-9) 

40 

. 1591549421 (-1) 

. 1 591549423 (-1) 

1 - 27 (-8) 

50 

. 1273239554 (-1) 

. 1273239538 (-1) 

7.04(-9) 

60 

.1061032940(-1) 

.1061032948 (-1) 

2.27(-8) 

70 

.0904568326(-2) 

.9094 5681 20 (-2) 

1 . 48 (-8) 

80 

.7957747219 (-2) 

.7857747102 (-2) 

1.47  (-8 ) 

90 

.7073552970(-2) 

.7073552977 (-2) 

9.42(-10) 

100 

.6366197815  (-2) 

.6366197677 (-2) 

2.17(-8) 

Normalized  for  E„  =  1. 

¥  /J#  O 

Note:  The  number  in  parenthesis  is  the  power  of  10  by  which  the  preceding 
number  must  be  multiplied. 
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TABLE  ( I V - 4 f ) 


Half-Volume  Check  of  Subroutine 
Versus  Eg  ( 68 ) 

Range  of  incidence  anglesi  0  5  0;$  ^0 
Distance  from  Edgei  o.ool  $  kf  $  *00 


e{ 

* 

ke 

Subroutine 

Value 

Eq  (58) 

Relative 

Error 

.  0 

0.001 

3-999999000 

4.0 

2.5(-7) 

-2.5 

0.001 

3.414213208 

3.414213562 

1.04(-7) 

*5.0 

0.001  . 

2.000000499 

2.0 

2.5(-7) 

r?.5 

0.001 

0.5857867902 

0.5357864376 

6.0(-7) 

90.0 

80.0 

6.7 57530467C-19) 

0.0 

6.76(-19) 

* 

Value  of  kp  where  the  maximum  relative  error  occurs. 

illative  error  here  is  defined  simply  as  the  difference 
between  the  Subroutine  and  Eq  (58)  because  one  cannot 
divide  by  zero. 
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Curves  from  top  to  bottom  are  for: 

e{  -  9c? j  MS0 
45:0 

&;s  2Z.S° 

6;  *  0.0° 


H - 1 - 1 - \ - 1 - 1 - 1 - 

.003  .005  .of 

k  p  (unitless) 


24.  Normalized  Power  Density  Vs.  k  ,TM 
Half-plane 


Vita 
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